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INTRODUCTION 


The present User's Guide applies to the three-dimensional viscous 
flow forward marching analysis as used for the calculation of the 
helicoptor tip vortex flow field. The basic analysis is described in 
Refs. 1 and 2 as well as Ref. 7 which present equations and results 
related to this program. The present User's Guide, which should be 
utilized in conjunction with Ref. 7 presents a discussion of the 
program flow and subroutines, as well as a list of input and output. 
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DESCRIPTION OF SUBROUTINES 


The present section describes the various subroutines used in the 
code. The main program is termed PEPSIG which controls the logic of the 
program. After setting initial conditions PEPSIG marches the solution 
downstream from the initial station. Subroutine PRIMRY in conjunction 
with COEFM and ADI solves the primary momentum equation. In this sequence 
COEFM writes the differential equation in finite difference form; i.e. 
loads the required storage locations to represent the finite difference 
form of the equation as a linearized algebraic equation at each grid 
point. This process yields a coupled set of linear equations, one 
equation at each grid point; which is then solved in ADI. Secondary 
flow is calculated via SECFLO; which calls COEFPH, ADI, COEFVS and 
ADI2X2. This set of subroutines serves to solve the continuity equation 
and the coupled stream function - vorticity equations. COEFPH forms 
the linearized continuity equation in terms of the unknown potential 
(see Ref. 7) at each grid point and ADI solves the set of equations 
through matrix inversion, COEFVS and ADI2X2 relate to the coupled stream 
function - vorticity equation. Of the remaining routines CROSEC is a 
routine structured to perform several unrelated objectives depending upon 
the argument list. These include Updating variables and turbulence quantities, 
computing integral quantities. Writing plot files, etc. RESTAR and RESTRG 
are restart subroutines. VISC and TURB set laminar and turbulent viscosity, 
respectively. A description of these and other routines follow. In addition, 
Fig. 1 shows the general flow through the program. 
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PRIMARY FORWARD 


PEPSIG 

Main Program 


Call INPUTS 
GEO 
IPROF 
CROSEC 


START, SET INITIAL 
FLOWFIELD 


START PROCEDURE 



March TO 
Next STEP 


FIG. 1 - OVERALL PROGRAM FLOW 
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ADDRES 


The first subscript in the Z array is the name of the variable stored. 
The second subscript is the grid point designator. The arguments of the call 
are the grid coordinates IY, IZ. The variables LN1A, LN2A and LN3A are 
passed through common from subroutine READZ, permitting subroutine address to 
operate both for the Z array incore option (INC0RE=1) or for the out of core 
option (INCORE=0). The output consists of the grid point designators used 
for Y direction derivatives KYI, KY2 , KY3 , Z direction derivatives KZ1, KZ2, 
KZ3, cross derivatives Kll, K12, K13, K21, K22, K23, K31, K32, K33 and source 
terms KPT. Note that shifting within ADDRES retains consistency of the 
output when used at boundary points. 


ADI 

Subroutine ADI solves 2-D scalar linear equations of the form: 

9(Aj>) = DU) + DU) + C 
3t y Z 

with boundary conditions of the form 

D y (*> + D zU) + C - 0 

where is a dependent variable, Dy and D z are Y direction and Z 
direction difference operators, and A and C are independent of <{>. A 
Douglas-Gunn scheme is used to solve the equations (Ref. 3). Parabolic 
equations are solved without iteration, elliptic equations are solved 
iteratively using a pseudo-time marching technique. 

The array C(N, IY, JZ) contains the initial conditions and the 
equations to be solved. The solution is returned in C( 1 , IY, JZ) for 
parabolic equations and in C(2, IY, JZ) for elliptic equations. 
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The argument list contains: 

M0DE1 

0 elliptic equation (A=l) 

-1 elliptic equation with no Dirichlet boundary 

conditions (A=l) 

1 parabolic equation 

2 paraboilc equation (A=l) 

IDUMP1 
0 
1 
2 

3 

4 

5 

INLY1, INLY2 
INLZ1, INLZ2 
MAXITS 

LIN 
ICON 

ADI2X2 


no output 

namelist and iterative information 
input C array plus IDUMP=1 

sweep 1 and sweep 2 coefficients plus IDUMP=2 

GAUSS coefficients plus IDUMP=3 

iterative information only 

index limits on IY 

index limits on JZ 

maximum number of iterations allowed 
for elliptic equations 

not used 

iterative convergence flag returned 
1 converged, 0 not converged 


Subroutine ADI2X2 solves two coupled linear two-dimensional equations. 
For parabolic equations the form is: 

(A f <|> f ) = D yf C<t> f . * s ) + D z f<*f» V " S fl *f " S f 2 + C f 

+ D yzf (4> f , <{> s ) 


( A s V = D ys (<f> f’ V + D zs (<l> f ’ V " S sl +f ' S s2 *s + C s 


3t 


+ D ( , 4> ) 
yzs y f s 


5 



For elliptic equations the form is: 

V*f’ V = D yf ( *f’ V + D zf <V * s > - s fi*f - S f 2^s + c f 

+ D yzf (< ^f * V 

V*f» * s > = D ys (< ^f ’ *s ) + D zs ( *f* V " S sl*f ” S s2*s + C s 

+ D y«s ( *f ♦.> 


where (j)f and 4> s are the dependant variables, D^-, Dy, D z are time, Y 
direction and Z direction operators. Dy Z is a cross derivative operation, 
and C is a constant. Subscripts f and s refer to first and second variables 
or equations. 

Boundary conditions are of the form: 

D y f C f » S ) + Dzf ( <J>f » <l>s ) + Cf = 0 


D ys($f » < l > s) + D zs(+f*+s) + C S “ C 

A Douglas -Gunn Linearized Block Implicit scheme is used to solve the 
equations (Ref. 3). Parabolic equations are solved without iteration, 
elliptic equations are solved iteratively using a pseudo-time marching 
technique. 

The array Q(N, IY, IZ) contains the equations to be solved. The array 
PHI(N, JY, JZ) contains the initial conditions and the solution. 

In the argument list IWRITE controls intermediate printout. Caution 
must be exercised with this variable since massive printout is produced. 
IWRITE 

0 iteration summary information only 

1 input equation dump plus IWRITE=0 

2 equations sent to block tridiagonal matrix inverter 
and solution returned plus IWRITE=1 

3 tridiagonal matrix inverter printout (subroutine 
MGS 2X2) plus IWRITE=2 

Argument ICONV is a flag set in ADI2X2 indicating convergence (IC0NV=5) for 
elliptic equations solved by pseudo-time iteration. Failure to converge is 
indicated by ICONV < 5. The maximum number of iterations allowed in 
iterating for convergence is MAX2X2 in common NAM2X2. Parameters MODEF and 
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MODES in the same common block indicate parabolic or elliptic equations and 
they are intended to be set to the same values for any problem, 1 for 
parabolic equations, 0 for elliptic equations with Dirichlet boundary 
conditions and -1 for elliptic equations with Neumann boundary conditions on 
all boundaries. 

ADI 3D (not used in the present problem) 

BULEEV (not used in the present problem) 

CL (not used in the present problem) 


COEFM 


This subroutine writes the streamwise momentum equation: 


0 = 


/ >u p 


<3u. 


h dx ■ 


‘ PV s 


/>u p w s ah 
~ ax’ 


au p 

ax, 


' P"s 


aup 

ax , 


P u p v s 


" T S7( p - tp «l 


dh 

ax. 


au r 


ah 


h aT — 


au r 


ah 


+ "H" axol ^ ax, ^ Up ax; 


ah au, 


+ -r- 


h ax, ax, 


/a U p / ah 

' ax 


/ \ 2 

ah au p /iUp / ah \ 

h ax 2 ax 2 h 2 l ax 2 1 


Calls to subroutine DXYZ apply the chain rule to express this equation in a 
body-fitted coordinate system used. Boundary conditions are either no-slip or 
symmetry as controlled by the indicies NS1, NS2, NS3, NSA for each cross- 
sectional boundary separately. 


7 



COEFPH 


The scalar potential equation: 



is written. Calls to subroutine DXYZ apply the chain rule to express this 
equation in body fitted coordinates. Boundary conditions set normal 
derivatives equal to zero on the wall and inner board and set <j>=0 at the 
center boundary. The resulting two-dimensional elliptic scalar equation 
in <)> is solved by a call to subroutine ADI. The arguments of the call are: 
MODE = -1 

IDUMP1 = ICOEF( 1,13) 

INLY1 = 1 
INLY2 = NEY 
INLZ1 = 1 
INLZ2 = NEY 
MAXITS = ICOEF( 1,3) 

LIN = 1 
ICONV = -1 


A warning message is printed if convergence has not been achieved. 

The irrotational components of the secondary flow velocity are computed 
from the scalar potential: 


3(J> 

H 

<ro 

1 

3<J> 

+ ^2 

3 4> 

= 8x l 

3x i 

3y l 

3x^ 

3y 2 

3 <f> 

= d l± 

3<t> 

+ ^2 

3<j) 

= 37 2 

8x 2 

3y l 

3x 2 

3y 2 


8 



COEFVS 


This subroutine writes the vector potential equation and the vorticity 
transport equation as a coupled system for solution by subroutine ADI2X2. 

Vector Potential 

dhi// \ d f I dh\j/ 
d x | J d x 2 l ph d x 2 

Vorticity Transport 
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Although zero normal derivatives of <f> in the scalar potential equation 
correspond to zero normal velocity, the tangential velocity, V fc , at solid 
boundaries is generally nonzero. In the coupled vorticity and vector-potential 
equations, both normal and tangential velocity components can be specified as 
boundary conditions since these equations are solved as a coupled system. 

By choosing (a) zero normal velocity and (b) -v t as the ^-contribution to the 
tangential velocity, the slip velocity v fc arising from the <f> calculation is 
cancelled, and the composite secondary flow velocity including both 4> and ip 
contributions will satisfy the no-slip condition exactly. Inner board and 
outer boundaries use boundary conditions of zero vorticity. The vector 
potential, ip , is set to zero on wall and inner board boundaries. At outer 
boundary, ^ is obtained by integrating inviscid vertical velocity. 

The equations are written in the array Q(N,IY, 1) for each line 
IZ=constant. Each line of equations is written onto logical unit NDR using 
binary write statements. These are to read in subroutine LOAD and used in 
SCALE, ADI 2X2 and MGS 2X2 to solve the equations. 

COPY 

The restart procedure saves all the common blocks in the COMDECKS COM 
and GCOM. In addition the entire Z array is saved as part of the restart. 
Whether the Z array is stored on disc or kept in core, subroutine COPY copies 
Z onto the restart unit or from the restart unit to the appropriate storage. 
Copy is called from subroutine RESTAR. 

CROSS (not used in the present problem) 

CROSEC 

Subroutine CROSEC is structured to perform several unrelated operations 
while surveying a cross section of the solution. This subroutine was written 
especially for the INC0RE=0 option to facilitate the use of auxiliary 
storage. This approach naturally inhibits vectorization. 
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Using the formal arguments I0PT1, I0PT2, I0PT3, IP0T4 subroutine CROSEC 
will perform up to four of the following functions at a single call: 


I0PT=1 

I0PT=2 

I0PT=3 


IOPT=4,5,6 


I0PT=7 


I0PT=8 


rotate N+l level quantities into N level storage 
update turbulence quantities 
compute integral properties 
unused 

load resultant secondary flow velocities and 
pressure at conclusion of computational step 

compute vorticity vector 


I0PT=9 load secondary flow velocities on starting iteration 

IOPT=10 write plot file and on M0DE=3 write geometry for 

potential flow 


DP (not used in the present problem) 
DFN‘ (not used in the present problem) 


DXYZ 


Subroutine DXYZ writes finite difference approximations in 
computational coordinates for terms of the form: 



( 1 ) 


where A and B are independent of <}> and Xj are reference coordinate 
directions. Transformation from reference coordinates to computational 
coordinates is performed within this subroutine. On option (LOAD) the 
resulting coefficients are loaded into appropriate storage for matrix 
inversion. Alternatively the term can be evaluated at the explicit (marching 
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step n) level. Setting derivative direction subscript i=0 produces terms of 
the form B 3 <f»/3Xj . Setting derivative direction subscript j=0 produces 
terms of the form A 3/3 ( B <J> ) . 

Marching direction derivatives are written as backward differences, 
transverse direction derivatives are written as central differences at the 
fully implicit (n+1) level or on option at the explicit (n) level. 

Transformation from reference coordinates, X, into computational 
coordinates, Y, is accomplished by the chain rule 3/3X£ = (3Yj/3Xi> 

9/8Yj with a summation on j. The elements of the Jacobian matrix are 
stored in the Z array. Difference weights in the transverse directions are 
stored in the AG array. 

The argument list contains: 


I,IJ 

computational indices of the grid point 

ND 

number of derivatives in summation in Jacobian 
transformation from reference coordinates to 
computational coordinates, ND=2 cross plane 
transformation only, ND=3 includes marching 
direction. 

ID , JD 

derivative direction subscripts (i,j) in Eq. 1. 

LOAD 

-i 

explicit evaluation of term for <j>=Z( IVAR, KPT) 

0 

difference approximation returned in AC array in 
common GRID 

1,2 

difference approximations loaded in Q array for 
ADI2X2. Equation is designated by LOAD (1 or 2), 
variable is designated by IVAR( 1 or 2) 

3 

difference approximations loaded into C array for 
solution by ADI 

IVAR 

(See LOAD) 

IDUMP 

0 

no printout 

* o 

printout of AC array 
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NPLUS 

0 Jacobian transformation evaluated at explicit (n) level 

1 Jacobian transformation evaluated at implicit (n+1) level 

Bet not used 


FRAME 


This subroutine, called from GEO, computes geometric quantities for the 
reference line based reference orthogonal coordinate system. Output 
includes : 


FFRAME( J, 1,1) 
FFRAME( J,2, 1) 
FFRAME( J,3, 1) 


normalized vector basis for 
the orthogonal reference 
coordinate system T, v, q 


B( 12, J, 4) 
B( 12, J ,5) 
B( 12, J,6) 


arc length derivative of the 
basis vectors 
3T , 3v , _9n 
as as as 


CNU.CETA 


components of centerline curvature 
K v » K r| 


A( 7 , 1) 
A( 7, 2) 
A(7,3) 


reference line arc length and 
derivatives 


B( 16, J,K) reference line locations in 

Cartesian coordinates and 
derivatives (k-1) 

Variables A and B appear in common GVAR, CNU and CETA are in common 
GEOS, and FFRAME is in common GEOM. In addition, B(12,J,K) contains the 
orthogonal basis vectors for the reference orthogonal coordinate system 
expressed in Cartesian reference coordinates (IFRAME-1) for use in computing 
coordinates for the potential flow or expressed in orthogonal reference 
coordinates (IFRAME=0) for use in the viscous flow. Thus, for IFRAME=0, 

B( 12, J,K) = 6 jk and for IFRAME=1 , B( 12, J,K)=FFRAME( J,K,1). The 

arc length derivatives of the basis vectors are also expressed in the same 

choice of vector bases depending on the IFRAME option. 
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Specification of the reference line is provided by either specification 

of the parameters S(t), <v(t), < n (t) on option IC0EF(1,9)=0 or 
specification of the location x(t), y(t), z(t) on option ICOEF( 1 ,9)=1. For 
either option, the following relations are used in integral (ICOEF( 1 ,9) a 0) or 
differential (ICOEF( 1 ,9)=1) form: 

3R = T 
3S 

AAA 

3T = K V + K n 

^ V n 


3v = 

as 


K T 

V 


dT) = - K 

Js 


T 

n 


where R(S) is the location of a point on the reference line and S is 
arclength. The following intermediate print is provided at each call to 
subroutine FRAME. 


WRITE(6,2) 

JX, [ (B( 16 , J,K) , J=1 , 3) , 

K= 1 . 3 ] 

WRITE(6,2) 

JX, [ ( A( 7 , J) , J= 1,3] 


WRITE(6, 1) 

FFRAME 


WRITE(6,3) 

JX,X,CNU,CETA 


WRITE( 6,1) 

[(B( 12, J ,K), J= 1 , 3 ) , K=4 

.6] 

1 FORMAT 

(' FRAME ' ,9F12.4) 


2 FORMAT 

(I6.9E11.4) 


3 FORMAT 

(' FRAME *, I12,3F12.4) 
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GAUSS 


Subroutine GAUSS solves a scalar matrix equation AX=B where A is a 
tridiagonal matrix and X and B are column matrices. Common GAUS contains the 
arrays, Cl, C2 and C3 which are the coefficients in the A matrix. The array 
C4 contains the entries of the B matrix. Gaussian elimination and back 
substitution are used to solve for X. The solution is returned in Cl. The 
indices II and IL indicate the first and last rows in the arrays to be used 
for inversion. 


GEO 

Subroutine GEO is the control subroutine for inputting geometry informa- 
tion and inviscid flow pressure and velocity. For the tip vortex generation 
problem, the geometry is created by the grid generation code of Frank Thames 
(Ref. 5). A separate subroutine GEOM reads and reprocesses the inputted 
grid system and computes Jacobian associated with grid systems. The 
reprocessed grid system and its Jacobian are stored in the coded file JOB30. 

The GEO subroutine will call subroutine TESTG which will read the grid point 
locations and its Jacobian from JOB30. At each computational grid point GEO 
loads the Z array with the transformation matrix from orthogonal reference 
coordinates to body-fitted coordinates, the determinant of the matrix, grid 
point locations in Cartesian reference coordinates. 

For viscous flow calculations (M0DE=1) subroutine GEO calls subroutine 
PFIELD to load the inviscid pressure field into the Z array. 

The first calf to subroutine GEO after a restart or at the initiation 
of a run reads namelist GEOM containing many of the parameters needed to 
specify the reference line and duct cross-section. Namelist input is explained 
in detail in the input section of this User’s Manual. 

INPUTS 


Subroutine INPUTS reads and prints namelist information. Details of 
the namelist input are given in the input section of the User’s Manual. 
Subroutine calls are made to VISC to initialize the laminar viscosity and to 
ROBTS to load the difference weights into the AG array and the computational 
coordinates into Y and YSAVE. 
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INTBL (not used in the present problem) 


1NVERS, INVERZ 

These subroutines invert 3x3 matrices. Some CDC computers require 
level 2 storage for some variables. INVERZ assumes the input matrix is in 
level 2 storage. 


IPROF 


This subroutine loads the initial conditions into the Z array. An 
inviscid velocity is computed from the input static pressure coefficient. 
Boundary layers are put on no-slip walls using subroutine VKPOHL for laminar 
flow and TPROF for turbulent flow. Laminar viscosity is loaded by a call to 

vise. 


LENGTH (not used in the present problem) 
LSFE (not used in the present problem) 


LOAD 


The equations for solution by the block coupled LBI solver, ADI2X2 , are 
written on logical unit LDR in subroutine COEFVS. Subroutine LOAD reads unit 
NDR and calls SCALE to perform time-step scaling to speed convergence of the 
solution. If resident memory storage limitation is a problem, the equations 
can be written on logical unit NDR2X2 in a form compatible with their use in 
subroutine ADI2X2. The controlling index is IN. The resident memory option 
is IN= 1 , and unit NDR2X2 is not used. Options IN=0 and IN=-1 both write unit 
NDR2X2 in the same format. Option IN=-1 requires less resident storage than 
IN=0, but requires more reading of unit NDR and consequently more computer 
time. 
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PEPSIG 


The code is based on the viscous primary /secondary flow analysis of 
Briley and McDonald (Ref, 1) as modified by Levy, Briley and McDonald 
(Ref. 2). The main program PEPSIG controls the logical flow or the code. 

The title card read at the top of the main program contains the index 
MODE in column 1. Successive cases may be read in, each with a new title 
card. If MODE=0 or the first four letters of the title are STOP, the 
program will terminate. 

M0DE=1 computes viscous flow using the subject viscous primary /secondary 
flow analysis. After initial geometry and fluid dynamics conditions are set or 
after restart conditions are read in, the main computational loop marches the 
solution downstream. Subroutine PRIMRY computes the primary flow by solving 
the streamwise momentum equation. Secondary flows, both rotational and 
irrotat ional , are computed by subroutine SECFLO. Post-processing and writing 
of a plot file are accomplished by calls to CROSEC. Printed output is 
provided through subroutine OUTPUT and saving restart file is controlled by 
subroutine RESTRT. The end of the loop is reached, thus finishing a marching 
step. Printed output, plot files and restart files are written at user 
designated stations and not necessarily produced at each marching step. 

MGS 2X2 


Subroutine MGS2X2 solves a matrix equation AX=B where A is a block 
triagonal matrix with block size 2X2. X and B are block column matrices with 
block size 1x2. The array Ql, equivalenced into common BLKMM, contains 
matrices A and B. Gaussian elimination, and back of substitution are used to 
solve for X. The solution is returned in matrix B. Arguments IL and IH 
indicate the first and last rows in the A and B matrices to be used for 
inversion. Argument IGDMP - 0 prints the input A and B matrices. 

NTRN 

Subroutine NTRN can be used to control transfers between fast core and 
large core memory on CDC machines. It can also be used to link with NTRAN 
utilities on Univac. 


17 



OUTPUT 


This subroutine controls the two-dimensional printout of variables 
from the Z array. The C array is used to transfer data from Z to the 
printout. The only argument of the call is: 

IWRITE 

1 printout of polar coordinates of each 

grid point, computational velocity components, 
streamwise vorticity, and pressure coefficient 

2 printout of velocity and vorticity for starting 
iteration. 

3 printout of auxilliary variables requested in 
ICOEF( 2 ,6 ) to ICOEF( 2 , 13) 

Velocity vectors are presented in computational coordinate directions. 
Additional printout of velocity in reference orthogonal coordinate directions 
can be turned on by setting ICOEF( 1 , 1 1 )=1 . 

The C array is loaded by a call to subroutine READ9 and two-dimensional 
output is provided by a call to subroutine 0UT2D. 

0UT2D 


This subroutine contains the write statements and formatting for 
two-dimensional printout of a single variable and a label. The label is 
passed through the argument list as MSG and the variable is selected from the 
C array by the argument IP. 
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PARSE (not used in the present problem) 


PFIELD 

A three-dimensional inviscid pressure field is input for use in the 
viscous primary/secondary flow analysis. The inviscid pressure field is 
supplied by Maskew [Ref. 6]. Subroutine PFIELD reads the inviscid pressure field 
input file, logical unit KDRUM and loads 0.5 times the inviscid pressure 
coefficient into the Z array. When no inviscid pressure field is being used 
IC0EF(1 , 6)=1 loads Cp=0 into the Z array. Standard operation sets IC0EF(1 ,6)=0. 

Since the inviscid flow can be run on a coarse mesh, streamwise smoothing 
of the inviscid flow pressure field is used to better represent the inviscid 
flowfield. The smoothed pressure coefficient at index JX is computed as 

CP(JX) = (CP(JX-l) + 2 *CP(JX) + CP(JX+l))/4 

PFLOW (not used in the present problem) 

POLY (not used in the present problem) 


PRIMRY 


The controlling subroutine for calculation of primary flow is PRIMRY. 
It calls COEFM to write the primary flow momentum equation, and calls ADI to 
solve the streamwise momentum equation. 

PROB (not used in the present problem) 

READZ, WRITEZ 

When the PEPSIG code is used on computers with very little resident 
memory it may be necessary to use extensive disc storage for many variables. 
On option INC0RE=0 the Z array containing the fluid dynamic and geometric 
quantities can be stored on logical unit 9. Subroutines READZ , WRITEZ and 
ADDRES work together to manipulate the Z array and keep track of storage 
addressing. READZ and WRITEZ call subroutine NTRN to interface with machine 
dependent utilities such as NTRAN or MOVLEV. 

For all operations subroutine READZ sets indices LN1A, LN2A and LN3A 
for use by subroutine ADDRES. 
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READ9 


For two-dimensional output of variables in a cross section, subroutine 
READ9 loads the C array for subsequent printout by OUTPUT and 0UT2D. The 
options associated with argument IWRITE are the same as those explained under 
subroutine OUTPUT. 


RESTAR, RESTRG 

A flexible restart capability is provided which allows the saving of 
the solution at designated stations and subsequent restart from those 
stations. The Z array, which may be stored on logical unit 9 (INCORE=0) and 
all the common blocks in common deck COM are written on the designated 
restart file JRSTOT or read in from the designated restart file JRSTIN. 
Subroutine RESTRG reads or writes the common blocks in common deck GCOM. 
Multiple restart files can be written on a single logical unit without 
multiple end of file marks. 


ROBTS 


Finite difference operators for first and second derivatives 
into the AG array in subroutine ROBTS. First derivative operators 


are loaded 
are: 



> 0 



and second derivative operators are: 



where h is the local computational grid interval. The computational grid is 
equally spaced going from 0 to 1. The computational coordinates are stored 
in the Y and YSAVE arrays by this subroutine. ROBTS is called by INPUTS at 
the start of a calculation. 
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SCALE 


A system of two coupled equations is solved by pseudo-time iteration 
using an LB I technique in subroutine ADI2X2. The time step taken at each 
grid point is scaled to help speed convergence. Four time scales are 
computed and the smallest time scale is used. The four time scales are: 

Tj twice the reciprocal of the geometric mean 

of the Y and Z direction operator diagonal terms 

T 2 16 * SIn2(tt/( 2 # NE)) divided by the magnitude of 
the source term 

T 3 4 * SIn2(it/( 2 *NE) ) divided by the magnitude of 

the difference between the off diagonal terms in 
the derivative operators 

T 4 a default value set to prevent a divide by 
zero ( 10 ^) 

SECFLO 


Subroutine SECFLO is the controlling subroutine for computation of 
secondary flows. Subroutine COEFPH is called to write and solve the scalar 
potential equation. COEFVS is called to write the equations for the coupled 
vorticity transport and vector potential equations. ADI2X2 is called to 
solve these coupled equations. The resulting solution is loaded into the Z. 
Special care is taken in computing the velocities at the polar coordinate 
singularity at the duct centerline. 


TESTG 

Subroutine TESTG is called by subroutine GEO to read the coordinates 
of grid point and its Jacobian from input file JOB30. 
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TPROF , VKPOHL 


For initial conditions a boundary layer profile is placed adjacent to 
a wall. Subroutine TPROF computes velocity profiles for a turbulent boundary 
layer based on a modified Coles profile. Subroutine VKPOHL computes velocity 
profiles for a laminar boundary layer based on a von Karman - Pohlhausen profile. 
These routines are called from IPROF at the start of a calculation. Density 
profiles are computed in IPROF based on velocity, pressure and the perfect 
gas law. 


TUBE (not used in the present problem) 
TURB 


TURB subroutine is called by CROSEC subroutine to update the eddy 
viscosity. The eddy viscosity computation is based on Claused Model, 
Eqs . (19 and 20) in the main report. 

UNIVAR 


Linear interpolation on a monatonic function is performed. The data is 
sampled to determine if the function is increasing or decreasing. The 
location in the data array of the previous data point is saved as a starting 
point for the new search. 


vise 


On an initial call from INPUTS with VISCOS=0, this subroutine sets the 
reference viscosity from reference density, velocity, length and Reynolds 
number. Subsequent calls from IPROF set laminar viscosity. 

WAVINT 


A weighted quadratic interpolation is performed using four node 
points. Quadratic equations are fit through the three points on the right 
and another quadratic is fit through the three points on the left. A 
weighted average of these two parabolas is provided. 
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DESCRIPTION OR INPUT 


TITLE CARD (I1,8A4) 

COLUMN 1 MODE = 1 viscous flow primary/secondary analysis 
MODE = 0 terminate run 


COLUMN 2-33 TITLE 


NAMELIST RESTAR 

IRSTIN is the station number to be read in for RESTART 

if IRSTIN = 0 not a RESTART case (DEFAULT IRSTIN = 0) 


IRSTOT is the interval for saving RESTART information 

If IRSTOT = 0 no RESTART info saved 


JRSTIN logical file name from which RESTART information 

is to be read (DEFAULT JRSTIN^ll) 


JRSTOT logical file name onto which RESTART information 

is to be written (DEFAULT JRST0T=11) 


NFILE is the sequence number in JRSTIN of the desired 

restart information (DEFAULT NFILE=0) 


NSAVED is the number of RESTART stations saved on JRSTOT. 

Must be initialized in inputs to the number of stations 
already written (and to be preserved) on JRSTOT. 
Nominally initialized NSAVED=0. 

NOTE that by setting JRSTOT = JRSTIN and NFILE = NSAVED 
one file can be used for both read and write without 
destroying the already saved information. 

INCORE Z Array uses no auxiliary storage (INC0RE=1), DEFAULT) 

Z Array stored on logical unit 9, retrieval and storage 
controlled by subroutines READZ and WRITEZ (INCORE=0) 
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NAMELIST FLUIDS 


X 

DX 

AP 

NFIRST 

NS 

NEY , NEZ 
IPLOT 

CMACH 

REY 

KTURB 

IBULEV 


marching parameter used in evaluation of polynomials 
P and PA in namelist GEOM. X is the streamwise 
computational coordinate. 

initial increment in X at station JX=NFIRST 


ratio of adjacent values of DX for geometric change 
in X. (default AP=1) 

station number of the parameter X. Note if this is 
not the same as the first station in a run subroutine 
INPUTS will compute the needed values of X and DX. 

station number of last station to be computed. 


number of grid points in circumferential and radial 
directions respectively. 

=1 plot file written as controlled by ICOEF (2,4) and 
ICOEF (2,5) or a controlled by IPSTA 
= 0 no plot file written 


reference Mach number 


reference Reynolds number based on UZERO, RZERO, YZERO 


= 0 laminar flow (DEFAULT) 

= 1 turbulent flow 

= 0 Buleev length scale not used (DEFAULT) 

= 1 Buleev length scale used in initial conditions 
to compute effective distance to wall only 

= 2 Buleev length scale used in calculation of 
turbulence mixing length 
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BLD(l) not used 


BLD(2) not used 

BLD(3) not used 

BLD(4) initial boundary layer thickness for IBULEV=0 option 

DELBUL value of Buleev length scale at edge of initial boundary 

layer for use with IBULEV=1 or 2 options 

ICDC indicator flag for disc writing in subroutine NTRN 

ARC centerline arc length at start of calculation 
(DEFAULT ARC=0) 

IWSTA station numbers printed output is to be produced 

(DIMENSION IWSTA (101)) activated for IC0EF(2,3) < 0 

IPSTA station numbers that plot file is to be written 

(DIMENSION IPSTA(lOl)) activated for ICOEF(2,5) < 0 

IC0EF( 1 ,N) 

N 

1 initial station print (0) 

2 momentum integral boundary layer calculation for boundary 
layer thickness used in turbulence model (1) 

3 PHI - ADI iteration limit (200) 

4 ADI3D iteration limit (200) 

5 max no. of pressure iterations in PRIMRY (5) 

6 0 - LOAD CPI, 1 - CPI = 0 (0) 

7 PFLOW INTER DUMP -1, FULL DUMP 1 or 3 (0) 

8y 

o _L printout from CROSEC (0) 

3x. 

J 
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9 

10 

11 

12 

13 

14 

15 

16 
18 
19 

ICOEF( 2 , N) 
N 
1 
2 
3 
A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

18 

20 


FRAME PA: 0 is CNU,CETA,S; 1 is Y,Z,X (0) 

OUTPUT PHI , PSI Cartesian velocities (0) 

OUTPUT Cartesian velocities (0) 

ADI2X2 dump (0) 

ADI dump for scalar potential equation (5) 

ADI dump for streamwise momentum equation (0) 
INPUT namelist dump (1) 

IPROF dump (0) 

Difference weights dump (0) 

PFIELD drum reading information (0) 


ADI2X2 iterations (100) 

PRINT INTERVAL (0) 

PRINT INTERVAL - NOTE: If .LT.O See IWSTA (1) 

PLOT INTERVAL (0) 

PLOT INTERVAL - NOTE: If .LT.O See IPSTA (1) 

AUXILIARY VARIABLE PRINTOUT 

AUXILIARY VARIABLE PRINTOUT 

AUXILIARY VARIABLE PRINTOUT 

AUXILIARY VARIABLE PRINTOUT 

AUXILIARY VARIABLE PRINTOUT 

AUXILIARY VARIABLE PRINTOUT 

AUXILIARY VARIABLE PRINTOUT 

AUXILIARY VARIABLE PRINTOUT 

N PLANE intermediate DUMP of velocity + vorticity (0) 
Geometry DUMP for MODE+3 (0) 


NOTE: Printout at station JX will be produced for ICOEF(2,3).GT.O If 

MOD(JX-ICOEF(2,2), IC0EF(2,3)) = 0 

The plot file will be written at station JX for ICOEF(2,5).GT.O 
If M0D( JX-ICOEF (2,4), ICOEF(2.5)) = 0 


Printout of auxiliary variables from the Z array are produced 
by filling ICOEF(2,6) - ICOEF(2,13) with the variable storage 
index from the following list 
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1 

u 

- velocity component at N 

2 

u 

- velocity component at N+l 

3 

V 

- velocity component at N 

4 

V 

- velocity component at N+l 

5 

w 

- velocity component at N 

6 

w 

- velocity component at N+l 

7 

p 

density at N 

8 

p 

density at N+l 

9 

ft 

streamwise vorticity component at N 

10 

ft 

streamwise vorticity component at N+l 

11 

c 

p I 

inviscid pressure coefficient at N 

12 

c 

p I 

inviscid pressure coefficient at N+l 

13 

J 

determinant of the Jacobian at N 

14 

J 

determinant of the Jacobian at N+l 

15 

y 

laminar viscosity at N 

16 

u 

laminar viscosity at N+l 

17 

ti T 

turbulence viscosity at N 

18 

y T 

turbulence viscosity at N+l 

19 

Y 

coordinate in v direction at N 

20 

Z 

coordinate in n direction at N 

21 

Y 

coordinate in v direction at N+l 

22 

Z 

coordinate in n direction at N+l 

23 

0 

dissipation function (2D:D) 

24 

£ 

turbulence mixing length 

25 


scalar potential 

26 

* 

vector potential 

27 

p 

pressure 

28 

- 

(not used) 

29 

V<J) 

scalar potential velocity component 

30 

w<}> 

scalar potential velocity component 

31 


vector potential velocity component 

32 


vector potential velocity component 
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33 


vorticity component in v direction 


34 vorticity component in n direction 

35 vorticity component in T direction 

36 Y reference Cartesian coordinate 

37 Z reference Cartesian coordinate 

38 X reference Cartesian coordinate 

39 - 47 components of Jacobian at N + 1 

48 - 56 components of Jacobian at N 

57 

58 

59 
60 
61 
62 


IC0EF(5,N) 

N 

5 modify pseudo-time increment for iterative 

solution elliptic equations 1/At = (l/At) REF 
2.**IC0EF (5,5) (DEFAULT ICOEF (5,5)=0) 


6 modify pseudo-time increment for iteration 

solution of coupled elliptic equations 
1/At = ( 1 / At ) REF 2.** ICOEF (5,6) 

(DEFAULT ICOEF(5,6)=0) 


modify pseudo-time Increment for iterative solution 
of 3-D elliptic equation 1/At = (l/AtWp 2.**IC0EF (5,7) 
(DEFAULT ICOEF(5,7)=0) 


time step selection in ADI3D is based on the magnitude of the 
finite difference operators. One time scale is based on a 
composite of the cross plane operators. Another time scale 
is based on either (ICOEF(5,8)=0) the streamwise operator or 
(ICOEF(5,8)=l) a composite of all three direction operators. 
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NAMELIST GEOM (not used in the present problem) 


File Input /Output 

The geometric information and inviscid pressure field needed in the tip 
vortex flow calculation for NACA0012 are stored in permanent files JOB30 and 
INVICI4. If restart calculation is desired, the restart information from 
previous calculations stored in local temporary file RESTEP should be saved 
and made permanent for restart calculation. Also if a plot is desired, the 
local temporary file TIPPL which contains plot information can be saved and 
made permanent. 

Sample Input Cards and Printed Output 

A sample case was run to illustrate the set up of inputs streams and 
typical printouts of the code. The given input card deck is for a new 
calculation. In the case of restarted calculation, proper parameters in 
namelist 1 RESTART ’ must be set as described in The User’s Guide. 

The code output first prints out the namelists, ’RESTART 1 , ’NAME’ and 
’FLUIDS’. This is followed by the computational coordinates of the grid 
points and their associated finite difference coefficients for the first and 
second derivatives in both directions. 

The next output item is the integral property across cross section, like 
the cross section area and mass flux. This is then followed by the printouts 
of the Cartesian coordinates of the grid points and of the flowfield variables 
at the initial starting step of the calculation. The enclosed output has been 
annotated to explain the printed variables, coordinates and iteration histories. 
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CONTROL CARDS TO RUN SAMPLE CASE 


IfUfll * T '•00 i C 67 0 0 0 * 

USER » 7 tSRAST:.., 

CHARGE ilOlEf / iLlvC • 

CATL I ST lLO = f ; i UN=795b3bC ) 
CATLISTtLO = F,UN = 7 f ‘5S3rC/ST = LPF} 
GET* INVICI A ) 

COPYE I ( IN V IC I A * REST I P3 ) 
CATALOG(kE STIP3*k»U»(\l) 

RETURN (INVICIA ) 


GE;T»0LCPL = PG5/UM=79553FC. 

UPDATE tP = Oi_C‘PL »L=A123A »F. 

TO VPS ( INPUT ,C6UD = C0,1PILL» REST I P3* UN = 7^5535 ♦ PW=SR ASTE , AC = 1 01 CG7) 


*/ FE.PSIG f-ODt. 

*1 EJL1102.1G 

XN=XPKtJPl> 

SJLTIP. 

USER(U=795535*PA=SRASTE,AC= 101867) 

RESOURCE (TL=A50, WS=153& ,LP=12, JCAT=MDBAT) 

FILES i PR I = * • 

ATTACH tCOMP ILE • 

ATTACH (RESTIP3) 

ATTACH (JOB30 > 

COPY( JOB30»GJOB) 

RETUR N ( JOB30 ) 

REQUEST(RESTIP/200D»RT=W) 

REQUEST(TAPE5»RT=W) 

COPY (INPUT »TAPE5) 

COPY (TAPl5*0UTPUT) 

£ REQUt ST* T IPPI /20Q0 .ftT = L > 

REQUf STi TAh L16/?000»RT=W) 

R EOUc S T ( T A PE 1 3/ 2 0 0 » RT = L ) 

REOUl S t t T APE 1<* / 2 0 CD » RT= w ) 

REOUE ST ( TAPE lb/20 00 , R 1 =W) 

F.lRTF: : AMl = OUTPU \ J 100 0 » I =COf*P IL.E ,H = B I WAR Y/2L0 , GPT = HL ) 
LOAD » f T f • A r> Y tL - OUTPUT t CN = G0 , 1 C U 0 , GRLF AL L = ) 

GO. 



NAMELIST INPUT FOR SAMPLE CASE 


u> 

N> 


1 S BtfiO TEST CASE -2 
P.RESTFT 
. IRSTH = G, 

IRST 0T=5 , 

JRSTOT=ll , 

JRST I N=1 1 » 

NSAVED=0, 

NFILE=0 

SENS 

S.FLUICS 

ITU8E=0, 

NEY=30,«EZ=33, 

I OVER DO , 

CHACH=Gloi *REY=ioOOOOO.OO*KTURB=l*ICC£Fat6)=i , ICOEF ( 1 *15) =1 

X=!UOOOotoX = lIoOOO ,AP=1.000G,WS=24,OELHUL=0.02 ? 

NS=5, 

NS1=0 ,NS2=0,NS3=0,NS*=0, 

IPUd=0, 

BLD C A )=0 • 070 , , 

ICOEF (2,15)=1,ICOEF(1,S>=0, ICOLF(l,16 ) -1. 

ICOEF (1,12> = C, , 

ICOEF (1,1 4) = 1, ICOEF (2, 17) = 1 » ICUEF ( 1 » I 3 ; = 1 * 

I PLOT -0 , 

ICOEF (3,1 ) = 10, 

ICOEF C3.2>=11 » 

ICOEF (3,3>=10, 

ICOEF (3,4)=14, 

ICOEF (3.5)=1, 

ICOEF (1,0 = 0. 

ICOEF (1 »lfe>=0 ♦ 

ICOEF (1,11)=1 » 

ICOEF (1,R)=0, 

ICOEF (2.16)=1 t 
ICOEF (1,12)=0, 

ICOEF (1 , 13 ) =1 » 

ICOEF Cl, 14>=1» 

ICOEF (1,10) = 0, 

ICOEF (2,1) = 100 , 

YS(1,1)=1. , 

YS(2,1)=30., 

YS (1 , 2) =1 • , 

YS (2 ,2)=33 • , 

XPP=1 ,2, 3, A , 

5, f 
21 . 

XPQ: 

0. 

0 

ol7814»0. 60 53 ,G»8273, 0.3450, 0.8635, 0.3377, 0.5 045, 

MFIRST=1 
SEND 
f.GEOh 
NGE0M=2, 

I TUB E =0 , 

ISYF=2 
&END 
STOP 



4455, 

7313, 



R*TMC TEST CASF 


PRINTOUT FROM SAMPLE CASE 


C 


& R F. S T k T 

IRSTIN 


u. 

IRSTGT 

- 

5, 

JR5TIN 

Z 

lit 

JRSTOI 

= 

lit ' 

NFILE 


0 N . 

NSAVED 

n 

0 • 

NS 



IfJCOKf 

&F.NO 


i 


NUN 

z 

1 % 

NIJ 

z 

2 ♦ 

NVN 

— 

3 . 

NV 

— 

‘1 

NVN 


5. . 

NV 

z 

c . 

NRMCti 

z 

7 . 

NR HO 

z 

‘ . 

NVORN 

Z 

0 . 

NV'OR 

z 

10 . 

NCI 5 1 N 

z 

1 i « 

NCR I 

Z 

let 

fJJN 

z 

1 : • 

NJ 

z 

14 i 

MUN 

z 

IF . 

MU 

z 

If . . 

M’JTN 

z 

17 , 

*UT 

z 

1-. 

■ J X Y Z N 

z 

I'm 

NXYZ 

— 

21* 

NDO 

z 

c' • 

NLtN 

z 

c f i t 

N 0 H I 

z 

?L». 

NPSI 

r 

'• . 

UPRES 

r 

? 7 * 

N V PH 

z 

2 V • 

MW PH 

z 

30. 

NVPS 

z 

31 * 

NVPS 

z 

o 2 ♦ 

NVOR1 

z 

7,3. 

NVCR2 

z 

34 , 

MV OP. 3 

z 

3b * 

NXYZA 

z 

35* 

ME 111 

z 

3S-* 

MENU 

z 

4 i . 

M Z VAR 

— 

74, 

MOV 

“ 

50. 

NDZ 

— 

50. 

NE1PF 

r 

10. 

112PF 

r 

10 . 

N I N 

z ■ 

5 1! . 

LZ 

~ 

It 5000* 

NSLAB • 

z 

3700* 

NREAO 

z 

0 . 

ME 

z 

50* 50 

N r F R K 
il’ND 

- 

5 0 



FLUIDS 

= .1E+01, 

X = .1E+01* 

S =5* 

P = .1E. + 01 » “ 

EY = 30, 

EZ = 33 , 

E1PF = 9* 

E2PF = 9, 

DUCT = 0, 

TUBE = 0, 

PLOT = 0, 

OVER = 0, 

MACH = .IE -01, 

EY = • 1 F + 0 7 , 

TURB = 1, 

LAU = 2, 

BULEV = 1, 

51 =0, 

52 =0, 

S!3 =0, 

SA =0, ' 

LD “ .0, .0, .0, . 7 f . - 0 1 , 

ELBUL = .2? -01, 


COEF = 0 

, 100, 

1C 

, c ^ 

c , c 

* 

n * ii, 

D* 0* 0, 

0, 0, 

0 

* 0 ♦ 

0 » t 

♦ 

Of 0 * 

O* 1, 0, 

0, 0, 

0, 0. 

0 

♦ It 

1, 0 

• 

0* 0, 

CDC = 1 

, 






FIRST = 1 

, 






S = . 

17 + 01, 

. 3: '. 

*U2t 

.13+01 

f 

* 2?>\+Q? 

PP = . 

1 • + 01 . 

r\ • 

r l. i. 

-01 * 

.37+01 

t 

. ^-ci . 

•14E+G2, 

. 15E+02, 

» 

1GV+02 

9 • 1 1 

p 

*02 t 

• 27 E + 02 , 

.2HE+02, 

♦ 

2 C : \ ♦ C 2 

» • 3 ■*. 

•■e 

0 2 « ,31 

.4E+C2, . 

0 , • 0 , 

a f; 

* • 0 • 

n 

% L- * 


* 0 • i+O « 

PQ = . 

1044E+G0 


,123 Or 

▼ CD* 

• 

1 4 0 4 ♦ 0 0 « 

.3312F+00, 

• 35o6 ■ 

+ U C* 

, . 3 • 

vs.xiu 

t 

. 4lf>l:\ •> 

'.6501E+00 i 

, 6777". 

*r c 

, 


« 

* 7313:. - 

•9045E+00 , 

» 0 , * 

cu 

.0, 

• 0 ! * 

0 

t • C A * 


END 


0. 0, 0, 0, in, C. C, 0, 0, 14, 3 , 0, 0, 0, 1, 

0, Ci« 0 • 0, 0. (j , 0, C , 0, 1, 0, 0, 0. G , 0, 0, 

0* 1» C • C » 0, 0 , 1 , . 0 , Ot 0, 0, 0, 0, 0, 0, 0, 


.S *01 • » ; ’ + U 1 » .77+01 , . °F + 01 , .--r+01, .1^+02, .111 + 0 

■ + 02* .1 •'••. + 02, .21+07, .21--+C2, .22C+C2, .23* + 0 2, .24*, 0 

•C.?. .32*1 < 02, .33'.+ 02, .34'+C2, .f=>G2. .307+02. .377 + 0 

*o* 7?n* • c , 

.1 GO, If + 00, .203 31+00 , .2276" +00* .25 If + 00, .2 

j , * fi ■*+ 1 *■ 0 G , .474'‘„+G0, i* ** 0 4 E ‘ '. + 0 0 , » o *> 1 * + 0 C , * tz ** 3 4 :. + Q 0 , 

1C. .7h7;-;*G0, . 7 1 4 - •» 0 0 , . >eHV/ + CG. .827?" + CC. .?4«* + CC* 

?• • C * * t) • ,0* .0, • 0 , *3, * C , *0, ,0, .0 
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COORD DUMP 
1 1.00000 
2 2.00000 

3 3.00000 

4 4.00000 

5 5.00000 


10 

11 

12 

13 

14 

15 

16 
17 
1 .-: 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 
1 
2 

3 

4 

5 

6 
7 
« 

9 

10 

11 

12 

13 

14 

15 

16 
17 
IF. 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 


6.00000 

7.00000 

6.00000 
9.00000 

10.00000 

11.00000 

12.00000 

13.00000 

14.00000 

15.00000 

16.00000 

17.00000 

18.00000 

19.00000 

20.00000 
21.00000 
22.00000 

23.00000 

24.00000 

25.00000 

26.00000 

27.00000 

28.00000 

2^.00000 

30.00000 
1.00000 
2.00000 

3.00000 

4.00000 

5.00000 

6.00000 

7.00000 

6.00000 

9.00000 

10.00000 
11.00000 
12.00000 

13.00000 

14.00000 

15.00000 

16.00000 

17.00000 

18.00000 

19.00000 

20.00000 
21.00000 
22.00000 

23.00000 

24.00000 

25.00000 

26.00000 

27.00000 

28.00000 

25.00000 

30.00000 
31 .0000 C 

32.00000 

33.00000 


- YF HYS » Y Ci> 

1 . 0 0 0 C 0 

r . o o o r (i 

3.00000 

4 . o 0 o r v 

5 . o o o r o 
r.oorno 

7 . o o c r o 

8 . o c o c c 

5.0000 0 

10*00 occ 
li.ooroc 
12.00000 

13.00000 
14 .00000- 

1 E . 0 C 0 c 0 

1 f . 0 0 0 0 0 

17.00000 

16.00000 

19.00000 

20.00000 
21.00000 
22. no COO 

23 .0 0 DCU 

24.00000 

25.00000 
2E .00000 
27 .00CO0 

28 .00000 

29.00000 
30 .oo ore 

1 . 0 0 0 0 c 

2.00000 

3.0000 0 

4 . 0 0 C 0 C 

6 . 0 0 C 0 0 
6 • 0 0 0 0 G 

7.00000 

8 . 0 0 0 C 0 
5.000CC 

10 .COCOC 
ll.OGCCG 

12 .ooonc 

13 .000 0 Cl 

14 .0000 0 

15.00000 

16.00000 

17.00000 

18 .0 0 0 0 0 

15.00000 

20.00000 
21.00000 
22.00000 
23.0 OCC 0 

24.00000 

25.00000 

26.00000 

27.00000 

26 .0 0 00 0 
26.00000 

30.00000 

3 1 * 0 0 G r - 0 

32 .00 0 0 0 
33 .00 C CO 


-.If OOOr +01 

-.foo no "+00 

-.5000 0^ + 1)0 

- . r n o n c ;- + o G 

-.‘7 0 0002 +00 

- . " 0 C 0 0 i. + 0 U 

- . * n o t.i o • . + o o 

-•5G00Gl+00 

- • r o o o o r + o u 

- . 5 0 0 0 C + 0 0 

-.*cocor+oc 

-»f 00 GO 1 :. + 00 
- .5 0 Ot'OL + GU 
-, r 0 OCOi+GO 
-.500003+00 
-.500000+00 

- . c 0 0 0 0 + 0 0 

- • 5 0 U 0 0 1 + 0 0 
-.* 0 0091 + 00 
-.'70000 i >00 
-.50000U+00 
- • * 0000 E + 00 
-.50000L+00 

- . v 0 t) l 1 0 l. + 0 0 
-.50000). + 00 
-. 50000 K + 00 
-.*00001+00 

- • 5 0 0 0 0 i\ +00 
-.5 0 00 0 f + 0 U 

• 50000 L’ + GO 
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The following three pages contain the values of the Cartesian vertical coordinate Y at the 
grid points (IY, IZ) in the first cross-sectional plane. Each column represents a given IZ location 
and each row represents a given IY location. For example, the ”Y" physical location for IY=10, 

IZ=4 is .05131. Following the print out for 'Y' , the Cartesian spanwise coordinate 'Z', the streamwise 
velocity 'u\ the vertical velocity 'v' , the spanwise velocity 'w\ the streamwise vorticity ’VOR-X' 
at grid point IY , IZ are printed out. 
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• 239 1E-0 1 
• 2444E-0 1 
•2511E-01 
•2596E-01 
•2704E-01 

• 284 2E- 0 1 

• 302 0E- 0 1 
•3250E-01 
•3550E-01 
•3947E-01 
•4481E-01 
•5214E-01 

• 6252E-0 1 
•7787E-01 


• 3157E-0 1 

• 3 159 E-0 1 

• 3161 E-0 1 
•3165E-0 1 
•3169E-0 1 
•3174E-0 1 
•3180E-0 1 
•3188E-0 1 

• 3198E-0 1 

• 321 1 E-0 1 
• 3227E-0 1 
• 3246E-0 1 
•3271E-0 1 
•3303E-0 1 
•3342E-01 

• 339 1 E-0 1 
•3453E-0 1 
•3532E-01 

• 363 1 E-0 1 
•3756E-01 

• 39 1 6 E-0 1 
•4120E-0 1 
•4383E-01 
•4723E-0 1 
• 5167E-0 1 
•5754E-01 
•6543E-0 1 
• 7627E-0 1 
•9162E-0 1 
•1143E+00 


LO 

00 


IZ= 21 22 23 24 25 26 27 28 29 30 

1 = *21005+02 « 22 00E + 02 • 2300E+02 •2400E+02 •2500E + 02 •26010E + 02 •270 0E+D2 *280QE+02 .2900E + 02 .3000E+02 


Y 

*1 000E+01 -•3958E-D1 -.45265-01 -.4 824E-01 -.49565-01 
•2000E+01 -.3960E-0 1 -.45305-01 -.4828E-01 -•4960E-01 

• 3000E + 01 -.3964E-01 -.4534E-D1 -.‘4833E-01 -.49665-01 
•4000r+01 -• 3968E-0 1 -•4539E-01 -.48395-01 -.4972F-01 
• 5000E+01 -.39735-01 -.45455-01 -.48465-01 -.49316-01 
. 60005 + 01 - • 3 9 6 0 E - D 1 -.4554 E-01 -.4S56E-01 -.49925-01 
•7000E+01 -.39885-01 -.45645-01 -.4863E-C1 -.50055-01 
.8000E+01 -.39995-01 -.45776-01 -;4883E-01 -.5022E-01 
•9000E+01 -.4012 E-01 -.45935-01 - « 49025-01 -.5043 r -01 

• 10Q0E + 02 - • 4 0 28E-0 1 -.46135-01 -.49255-01 -.5069E-01 
•11005+02 -.404 9 E-01 -.4638E-01 -.49555-01 -.5102E-O1 
« 1 20 0 E+ 0 2 - . 40755-D 1 -.4670E-01 -.49925-01 -.51445-01 
.13QCE+02 -.41085-01 -.47105-01 -.50335-01 -.51965-01 
•1400 E +02 -•41485-01 -.47595-01 -.5096E-01 -.52615-01 
•15005+02 -.42005-01 -.46225-01 -.51705-01 -.53435-01 

• 1 600E + 02 4 264 E-0 1 -.49015-01 -.5262E-01 -.54465-01 
•1700E+02 -.4345E-01 -.5000E-01 -.537HE-01 -.5576E-01 
• 1PO0E+O2 -.4448E-01 -.51265-01 -.5524E-01 -*57405-01 

• 1900E + 02 -.45785-01 -.5284E-01 -. , 5709E-01 -.59485-01 
*20005+02 -.47425-01 -.54855-01 -.55445-01 -.62115-01 
.21005+02 -.49515-01 -.57415-01 -.62425-01 -.65465-01 
.22005 + 02 - .52.1 6 E-0 1 -.60675-01 -.6624E-01 -.69735-01 
•23005+02 -.55625-01 -.64675-01 -.71145-01 -.75235-01 
*24005+02 -.60075-01 -.70305-01 -.77495-01 -.82345-01 

• 2500 E + 02 -.65875-01 -.77405-01 -. 85785-01 -.91635-01 
.26005+02 -.73565-01 -.86795-01 -.96745-01 -.10395+00 
.27005+02 -.83885-01 -.99415-01 -.11155+00 -.12045+00 
♦28005+02 -.98065-01 -.11675+00 -.12175+00 -.14315+00 

• 29005*r02 -.11815 + 00 -.1.4156+00 -. 16046 + 00 -.1753000 
.30005+02 -.14785+00 -.17756+00 -.20275+00 -.22276+00 


.49935-01 -.49995-01 -.5000E“01 -.5D0DE-01 -*5000E-01 
.49985-01 -.50045-01 -.5005E-01 -.5005E-01 -.5005E-01 
.50045-01 -.5010001 -.5011E-01 -.5O11E-01 -.5011E-01 
.50115-01 -;5018001 - • 50 19 E- 0 l -.5O19E-01 -*5019E-01 
.50206-01 -.50285-01 -.50295-01 -.50295-01 -.5029E-01 
•5032E-01 -.50405-01 -.5041E-01 -.5041E-01 -.5041E-01 
.50465-01 -iBOSS^-Ol -.50565-01 -.5057E-01 -.5057E-01 
.50645-01 - ; 50745-0 1 -.50765-01 -*50765-01 -.50765-01 
.50875-01 -.50985-01 -.510CE-01 -.51005-01 -.51005-01 
.51165-01 -.51275-01 -.51305-01 -.51315-01 -.513-15-01 
.51515-01 *-.51655-01 -.5168E-01 -.51695-01 -.5169E-01 
.51965-01 -452125-01 -.5216E-01 -.52175-01 -.52175-01 
.52525-01 -.52705-01 -.52765-01 -.5277E-01 -.52775-01 
.53235-01 - i 53445 -0 1 -.53515-01 -.5353E-01 -.5353E-01 
.54125-01 - i 54375-0 1 -.54455-01 -.54485-01 -.54485-01 
•5523E-01 -.55535-01 -.5564E-01 -.55675-01 -.5567E-01 
.56645-01 -.57005-01 -.57145-01 -.57185-01 -.57185-01 
.58415-01 -458865-01 -.5903E-01 -.5908E-01 -.59085-01 
.60655-01 -.61205-01 -.61426-01 -i6148E-01 -.614-8E-01 
*63495-01 -464175-01 -.6444E-01 -.6452E-01 -.6452E-01 
.67115-01 -.67955-01 -.68295-01 -.63395-01 -.6839E-01 
.71735-01 -.72785-01 -.73225-01 -.73345-01 -.73345-01 
.77675-01 -.78995-01 -.79545-01 -.79705-01 -.7970E-01 
.85365-01 -487035-01 -.87735-01 -.-87945-01 -.87935-01 
.95405-01 -.97535-01 -.98435-01 -.98695-01 -.9869E-01 
.10875+00 -.11145+00 -.11265+00 -.1129E+00 -.11295+00 
.12655+00 -il3015+C0 -.13165+00 -.13205+00 -.13205+00 
.15115+00 -.1557-+00 -.15775+00 -.1583E+00 -.15835+00 
. 18566 + 00 -.1920-+00 -.19475 + 00 -.195.55 +CO -.19555+00 
.23715+00 -.24575+00 -.24935+00 -.25045+00 -.25045+00 


-.5000E-01 
-.50055-01 
-.50115-0 1 
-.50195-01 
-.50295-01 
-•5041 E-01 
-.5057E-01 
-.50765-01 
-.51005-01 
-.51315-01 
-.51695-01 
-.52175-01 
-.52775-01 
- • 5352E -0 1 
-.54475-01 
-.55665-01 
- . 5717E-0 1 
-.59065-01 
-.61465-01 
-.64505-01 
-.68375-01 
-.73315-01 
-.79665-01 
-.87885-01 
-.98615-01 
-.11285 + 0 C 
-.13195+00 
-.15815+00 
-.19535+00 
-.2501E + 0 C 



IZ= *31 32 33 

Z = #31002+02 . 32 0 0E + 02 .33 ODE* 02 

1 # 1000 E + 0,1 -.49992-01 -.45992-01 -.5000E-01 

2 #20002+01 - #50052-0 1 -.50052-01 -.50052-01 

3 • 3 00 0 E + 0 1 -.501U-G1 -.50112-01 -.’50112-01 

4 .4 000E+01 -.50158-01 -.50198-01 -.50192-01 

5 *50002+01 - • 5 0298-D 1 -.50292-01 -.50292-01 

6 .6000E+01 -.50412-01 -.5041E-01 -..5041E-01 

7 •7000E+D1 -.50562-01 -#5056 E- 01 -.50568-01 

8 * 8 00 OE +0 1 -.50762-01 -.5076E-01 -.5076E-01 

9 .9000E+01 -.510QE-D1 -.51008-01 -.5100E-01 

10 .10008+02 -.51318-01 -.5131E-01 -.51312-01 

11 .1100E+02 -.51692-01 -.5169E-01 -.51692-01 

12 . 1200-E + 02 -.52171-0 1 -.52172-01 -.5217E-01 

13 » 1300 8 + 02 - • 5277 2-01 -.52772-01 -.52772-01 

14 .14008+02 -« 5352E -0 1 -.53528-01 -.53522-01 

15 .15008+02 -.54478-01 -.54478-01 -;5447E-01 

16 ,16008+02 -,55668-01 -.5566E-01 -.55662-01 

17 • 17008+02 -,57161-01 -.57168-01 -.5716E-01 

18 .18008+02 -.59058-01 -.5505E-01 -.5906E-01 

19 .1900E+U2 -.61452-01 -.61458-01 -.6145E-01 

20 ,20008 + 02 -*64498-01 -.64492-01 -. 6449E-01 

21 .21008+02 68352-01 -♦66352-01 -.'68362-01 

22 .2200E+02 -.73252-01 -.73298-01 -.7329E-01 

23 .2300E+02 -.79638-01 -.79632-01 -.7964E-01 

24 .2400E+Q2 -.87852-01 -.8785E-Q1 -.8786E-01 

25 .25002+02 -*98588-01 -.9858E-01 -.9 859E-01 

26 .2600E+02 -.11281+00 -.1128E+00 -^11282+00 

27 . 2700E +02 -.13182-00 -.13182+00 -.1319E + 00 

28 .2800E + 02 -.15802*00 -.158QE+GG -.1581E + 00 

29 • 2900 E + 02 -.l q 518«DG -.19511 + 00 -.19528+00 

30 .3000E+02 -.24958+00 -.24592+00 -.25002+00 


STATION 


» 1 000 E + 01 .2GG08+01 • 30008+01 


4 5 

.4000 E +01 .50008+01 


. 60 OJD £ ♦ 0 1 *70 0 0801 .8000E + 01 .9000E + 01 


.10008+01 
.20008+01 
.30002+01 
.40008+01 
.50005+01 
.60005+01 
•70002+01 
.acoor+oi 
•90002+01 
.10001+02 
. 1100E+02 
•120CE+02 
•1200E+C2 
.14002+02 

• 1 50 0 r. + 0 2 
•16008+02 
. 17002+02 
.ifoor+02 
. 19002+02 
♦2000 E +02 
.2 1008 + 02 
.22001+0? 
.23001 + 02 
•2400E+C2 
•250CE+C2 
.26002+02 

• 2 7 0 0 8 * 0 ? 

• 2*.C0r + 02 

• 2 c . 0 0 >02 

♦30001+02 


. 25008-00 
.25001+00 
*2500 r +DO 
.25001+00 
.25001+00 
.25DOr+DO 
.25005+00 
•2500 L +00 

• 25 0 OE + O 0 

• 2 5*0 0 E ♦ 0 0 

• 25 o or +00 

.2500? +00 

• 250 CE + 0 0 
•2500 l +00 

• 250 OE+O 0 
•25008-00 
#25008*00 
.25002+00 
.2500^ +00 
.25002+00 
.25002-00 

- • 25 0 08 + 0 0 
- .250 CE + 0 0 
-• 25U0E + 0 0 
-.25002+00 

- * 2 5 C 0' -00 
-.25002+00 
-.2 r :0C‘ +00 
-.2500*+ Of 

- • 2 5 D 0:>C 0 


•21058+0 G 
•2105E+00 
« ? 1058+00 

• 2 105 E + 00 
•21G5E+00 

• 2 105E + 0G 

• 2 1 05E + 0 0 
•21058+00 
. 2 105E + 0 0 
•21U5E+00 
.21 0 5E + 0 0 

• 7 1U5E+00 

105E + 00 
*2 105 + 00 
•?105l+00 
•2105E+00 

• 210 58 + 0 0 
.2 1C5E+GG 
.? 10 52+00 

• ? 1 C 5 r. + 0 0 
.21042+00 

• 2 1 0 A E + 0 0 
.21042+00 

• ? 1 0 4 E+ 0 G 
.210^E+00 

• 2 1032+00 
.21032+00 

1022+0 c 

- . r 1 U 1 :, + 0 0 
-•210CL+00 


. 1 7542 + 00 
.17542+00 
. 1 7542+00 
.17548+00 
.17542+00 
.17542+00 
•17542+00 
•1754E+00 
. 1 754E+00 
.17542+00 
. 1 7542+00 
. 175*2*00 
; 1 7 5 3 2 + 0 0 
.17538+00 
. 1753E+00 
. 1753E + 00 
« 1 752E + 0 0 
. 1 7 t= 22 + U 0 
. 17512*00 
.'-1 7502 + 00 
. I74 q 8+0G 
,1 748E+00 
, 1 7462+00 
17442+00 
-.17412+00 
-.17372+00 
17322+GO 
1 7?2;> C 0 
-.171 5 2 + 0 U 
-. 1700L + 00 


•14432+ 00 
•14428+00 
•14422+00 

• 14428+00 
.14428+00 
#14422+00 
.14422+00 
.14422+00 
. 1442E+ 00 

• 1442E+00 
. 14415+CO 
. 1441E+00 
•14412+00 
•14402+00 
•1439&+00 
•1438E+00 
•1437E+00 
« 14 362+00 
•14342+CO 
.14322+00 
♦14298+00 
•14268+00 
.14218+00 

-*14158+00 
-.14058+00 
13952+ CO 
-♦13^42+00 
-.13 i. 5 E +00 
-. 1339C+00 
-.13008+00 


•#11652+00 
-.11658+00 
-.11652+00 
11658+00 
-.11648+00 
-.11648+00 
-.11642+00 
-•11648+00 
-.11638+00 
-. 11638+00 
-•1163E+00 
-.11625+00 
-.11618+00 
-.11602+00 
-.115SE+00 
-.11572+00 
-.11552+00 
-.11538* 00 
-.11502+00 
-*1146 r +00 
-•11408+ GO 
-.11342+00 
-.11265+00 
-•13 158+00 
-.11005+00 
-.10328+00 
-.10568+00 
-.10225 *00 
-.97245-01 
-. P997E-01 


«9l69~-01 
•91672-01 
.9166^-01 
•S164T-01 
•91622-01 
. 91592-01 
J9156E-01 
. 9 1528-0 1 
.91472-01 
. 9 1 4 1 ^ - 0 1 
.91332-01 
.91232-01 
•91102-01 
.90952-01 

• 9 07 58-0 1 
.90502-01 
.90198-01 

• B97 98-0 1 
.89292-01 
. 86662-01 
•87852-01 

4 86822-01 

.85502-01 
-.83782-01 
- « 8 1558- C 1 
785.9~-Q 1 
-.74612-01 
-.6*114- -01 
-4614 02-0 1 
-.49978-01 


-•69462-01 
-.69452-01 
-•69438-0 1 
-.65418-01 
-.69382-01 
- • 6934 E-0 L 
-.6929E-01 
- • 69 24 E-0 1. 
-.69162-01 
-.65072-01 
-.68968-01 
-.68822-01 
-.6864T-01 

-.68412-01 
-.68132-01 
‘•.6778 2-01 
-.67332-01 
-.6677 2-01 
-.66052-01 
-.65152-01 
-.64002-01 
-.62532-01 
-.60642-01 
-.58202-01 
-.55008-01 
- • 50 76 8-0 l 
-.45102-01 
-.37302-01 
-.26268-0 i 
-.^9508-0 2 


.49488-01 
.49468 -01 
.49438-01 
.49402-01 
.49368 -01 
.49322-01 
.49252-01 
.49188-01 
.49082-01 
.48968-01 
.48818-01 
*48622-01 
•4838^-01 
.48082-01 
,47708-01 
.47232-01 
.46632-01 
.45802-01 
.44932-01 
.43728 -01 
.42182-01 
• 4D 22 E-0 1 
.37702-01 
.34442-01 
.30 178-01 
-.24532-01 
-.16958-01 
-.6^378-0? 
.82062-02 
•2 q 99 8-01 


.31422-01 
.31392-01 
.31368-01 
.31328-01 
.31272-01 
.31218-01 
.31132-01 
.31038-01 
.30912-01 
.30768-01 
.30578-01 
.30332-01 
.30032-01 
• 2965E* -0 1 
.29178-01 
.28572-01 
.27828-01 
•26862-01 
.25662-01 
.24132-01 
-.22192-01 
-.19718-01 
-.16528-01 
-.12392-01 
-#69942-02 
.14098-03 
.97308-02 
.22908-01 
.41558-01 
.69118-01 


-.15002-01 

-.14978-01 

-.14932-01 

-‘.14888-01 

-.14828-01 

-.14752-01 

-.14658-01 

-.14542-01 

-.14398-01 

-.14202-01 

-.13972-01 

-.13688-01 

-.13322-01 

-.12862-01 

-.12288-01 

-.11562-01 

-•1065E-01 

-.95008-02 

-.80462-02 

-.62028-02 

-.38562-02 

B5868-03 

#29932-02 

.79792-02 

.14492-01 

.23112-01 

.34692-01 

.50602-01 

.73112-01 

.10648+00 




IZ= 

11 

12 

13 

14 


Z = 

« l'l QO £■*• 02 

• 12 0 0C+ 02 

•1300E+02 

r •1400E+02 

IY 

1 

Y 

•.lOOOE + Ol 

•3420E-D5 

• 1 36 3E-0 1 

♦ 2549E-01 

• 35 0 1 E- 0 1 

2 

•2G00F+01 

.3891E-04 

• 1 36-7E-0 1 

* 2 553E-0 1 

.35055-01 

3 

• 3 000 E + 0 1 

.8332004 

. 1 372E-0 1 

•2558E-01 

• 35 1 IE-01 

4 

• 4 000 E + 0 1 

• 1389E-0 3 

• 1 378 0 0 1 

• 2 565E-0 1 

•3518E-01 

5 

•5000E+01 

• 2085003 

• 1 386E-0 1 

• 2574E-01 

•3527E-01 

- 6 

• 6 OOOE + Ol 

•2957E-03 

•139&E-01 

•2584E-01 

•3539E-01 

7 

.7D00E + 0:i 

• 4 0 4800 3 

• 1 406E-01 

• 2 59GE-0 1 

. 3553E-0 1 

8 

•8 OOOE+Ol 

•541SE-03 

. 1 423E-0 1 

•261SE-C1 

•3571E-01 

9 

• 9 00 0 E + 0 1 

• 7 127E-D 3 

•1443E-0X 

•263 6 E-01 

•3593E-01 

10 

•1000E+02 

•9272E-03 

•1467 E-01 

• 2 662E- 0 1 

•3621E-01 

11 

•1100E+02 

•1196E-02 

• 1 45.7E-0 1 

. 2 69 5E-0 1 

• 3657E-0 1 

12 ’ 

•1200E+02 

• 1 533E-D 2 

•1535E-01 

•2736E-01 

•3701E-01 

13 

• 1 300 E.+ 02 

.1956E-02 

•1582E-01 

•2788E-01 

.37565-01 

14 

•1400E+02 

•2488E-02 

•1642E-01 

*2853 E-01 

• 38 26E - 0 1 

15 

• 1 500E + 02 

•3155E-02 

•1717E-01 

•2535E-01 

•39 13E-01 

16 

. If aoe+02 

• 39955-0 2 

.iniE-oi 

.3 03 RE- 01 

•4023E-01 

17 

• 1 70 0 E + 02 

.5053002 

- 1 S30E-0 1 

• 3 167E-0 1 

•4161 E” 01 

18 

•1800E+02 

• 6387E-D 2 

.? 07 9 E-0 1 

• 3 331 E-0 1 

• 4336E-0 1 

19 

•1900E+02 

i 80 75E-D 2 

•2269E-01 

•3538E-01 

•45575-01 

20 

• 2 000E + 02 

•1022C-0 1 

•2509E-01 

. 3800 E ”01 

•46375-01 

21 

• 2 100E + 02 

•1294 E-01 

. 2 o 1 4 E - 0 1 

•4133E-01 

•51535-01 

22 

♦ 2'? 0 0 E + 0 2 

•1642E-01 

•2205E-01 

•4E60E-01 

.56485-01 

23 

•2300E+02 

•20S9E-01 

♦3706 E-01 

•510 8 E-01 

•6233E-01 

24 

• 2400E + 02 

•2668E-01 

.4 356 E-0 1 

•5817E-01 

• 6991E-01 

25 

• 2 SClOE + 02 

•3424 E-01 

•5204E-01 

•G743E-01 

.79805-01 

26 

• 2600E+02 

•4424E-D1 

•6326E-01 

.75686-01 

•92885-01 

27 

♦2700E+02 

•5768E-01 

. 7 634 E-0 1 

.9615E-01 

.11 05E+00 

28 

•2EOOE+02 

•7614E-01 

.9906E-01 

.1 188E + 00 

•1346E+00 

29 

•2900E+02 

•1023E+00 

.1234E+00 

• 1508E+00 

•1688E+00 

30 

• 3000E+02 

• 1409E + DCI 

• 1717t>0C* 

« 1 581E+0 0 

• 21‘94 £ + 0 0 


IZ= 

21 

22 

. 23 

24 


Z = 

• 21 0 0 E + 0 2 

• 22 0 0 E^O 2 

• 230 0E + 02 

•2400E+02 

IY 

1 

Y 

• 1 00 0 E + 0 1 

•2549E-0 1 

• 1 363E-0 1 

•3420E-05 

- .1500E-01 

• 2 

• 200 0E + 0 1 

•2553E-01 

•1367E-01 

• 3 89 1 E- 04 

-•1497 E-01 

3 

• 300 0 E + 0 1 

• 2558E-0 1 

• 1 37 2 £- 0 1 

•8332E-04 

-• 1493E- 01 

4 

• 4 OOOE + Ol 

♦2565E-01 

• 1 37 8 E - 0 1 

• 1 389E-0 3 

-.1488E-01 

5 

• 5 00 0 E + 0 1 

•2 574 E-01 

•1386E-01 

•2085E-03 

- • 1 4.82E- 0 1 

6' 

•6000E+01 

•2584E-01 

•1396E-01 - 

•2557E-03 

-•1475 E-01 

7 

•7000E+01 

.25985-01 

•1408E-01 

« 4 0 4 8 E - 0 3 

- • 14 65E- 0 1 

8* 

•8000E+01 

• 2 61 5 E-0 1 

• 1 423E-0 1 

•5415E-03 

-•14 54 E-01 

9 " 

’ *9 OOOE+Ol 

4 2636 E-0 1 

•1443E-01 

•7127E-03 

-•1439E-01 

10." 

* 1 000 5"+ 02 

• 2662E-D 1 

•1467E-01 

•9272E-03 

- • 14 20E-01 

11 

•1100E+02 

•2695 E -01 

• 1 497E-01 

. 1 196E-0 2 

- • 1397E-0 1 

12 

• 120QE+02 

• 2736E-0 1 

• 1 535E-0 1 

• 1 53 3 E-0 ? 

-.1368E-01 

13 

• 1300E+02 

• 2788E-D.1 

•1582E-01 

•1956E-02 
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5 *5000 E + 01 .0 

6 .6000E+01 .0 

7 .7OO0E+O1 .0 

a .aoooE+oi *o 
9 .9000E+01 *0 

10 .1000E+02 .0 
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.2^00 f + 02 
•2900 r +02 
•3000E+C2 


• 35L,75“lb 
. 9 3 6 7 e - G 9 
.4 6665— 0 8 
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1 • 1 0006+01 .6033E-04 .14716-03 .43716-03 .4061E-03 

2 .2000^+01 -.2236E-D3 -.45746-03 -.83866-03 -.11786-02 

3 .30006+01 -.46396-03 -.1122E-02 -.1881E-02 -.28226-02 

4 .4000L+01 -.79586-03 -.17946-02 -.31736-02 -.47796-02 

5 .50006+01 -.11426-02 28136-02 -.48296-02 -.73526-02 

6 .60006+01 -.16176-02 -.29516-02 -.72456-02 -.1060E-01 

7 . 7 00 0 E + 01 -.2364E-02 -.5317E-02 -.98006-02 -.1442E-01 

8 .80006+01 -.317OE-02 -.7207E-02 -.12876-01 -.1903E-01 

9 . 900 0 E + 01 -.4054E-02 -.9707E-02 -.1724E-01 -.2482E-01 
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11 . 1 100E+02 -.71726-02 -.16236-01 -.25176-01 -.4289E-01 
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23 

•2300E+02 

-.2187E+00 * 

•4099E+00 - 
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5 . 5 000 E + 0 1 -.4829E-02 -.2813E-02 -. 1 142E-02 -.3494E-03 
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7 . 7 P00E + 01 -.9800E-02 -.5317E-02 -.2364E-02 -.6942E-03 

8 .GOOOE+Ol -.1287E-01 -.7207E-02 -.3170E-02 -.93976-03 

9 .9000E+01 -.1724E-01 -.9707E-02 -.4054E-02 -.1245E-02 
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13 . 1300 E+0 2 4843E-0 1 -.2737E-01 -.120eE-01 -.3775E-02 

14 .1400E + 02 - . 62 15E-0 1 -.3525E-01 -.1574E-01 -.4954E-02 

15 .15006*02 -.7960E-0 1 -.4545E-01 -.2054E-01 -.6637E-02 

16 .1600E+02 -.10166+00 -.6905E-01 -.26876-01 -.8805E-02 
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*11026 + 00 -. 1239E + 0 0 -.1296E + 0D -.12396+00 -.11026+00 -.8961E-01 
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.46316+00 -451276+00 -.5309E+03 -.51276+00 -.4631E+00 -.3836E+00 
.58376+00 -.65126*00 -.67286+00 -.6512E+00 -.5887E+C0 -.4912E+00 
• 7522E+00 -.82756+00 -.85416+00 -.82756+00 -.75226+00 -.6327E+00 
•9636E+00- -.10566+01 -.1089E+01 -.1056E+01 -.9636E+00 -.8190E+00 
. 1240E+ 01 -.13546 + 01 -.13936 + 01 -.13546*01 -.12.4-0E + 01 -.1064E + 01 
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.34166-05 *42016-06 
.15256-04 -*20436-05 
• 3333E-04 - . 49846-05 
.53266-04 -.95036-05 

• 80.58E-04 -.13045-04 
.11276-03 -.17256-04 
.15896-03 -418576-04 
.22656-03 -.36046-04 
.28826-03 -.43736-04 
•3871E-03 -.68436-04 
.49606-03 -.80166-04 
•6712E-03 -.10476-03 
.9248E-03 -’.13666-03 
.1201 E-02 -.20915-03 
. 1 59 4 E- 02 - *27846-03 
.22126-02 -.35176-03 
.2933E-02 -456576-03 
•4074E-02 -.77626-03 

• 571 IE-02 -. 10905-02 
.79486-02 -.16365-02 
. 1153E- 01 -.23415-02 
.16576-01 -.35565-02 
.24036-01 -.53255-02 
.35415-01 -.81665-02 
. 520 IE- 0 1 -. 12555-01 
.12446-01 • 48225-0 1 
.92956-01 -436395-01 
.0 .0 

.0 .0 

•7615E-14 441105-14 


.18246-07 -.6303E-08 
•13596-09 .32075-08 

.5162E-0 9 .3430E-06 

.5450 E-OS .1364E-05 
.14756-09 .19496-05 

.9877 E-0 9 -.62365-06 
.18155-05 -.12125-05 
.85575-09 -.22596-05 
. 5924 E-0 5 -.4053E-05 
.52 18 E-0 5 .3011E-05 

.61116-05 41291E-05 

.1171 E-0 4 -.2730 E -05 
. 3690 E-0 5 -.10686-04 
. 5280 E-0 5 .20596-05 

.14886-04 . 1868E -04 

•17606-04 -.43136-05 
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.28956-04 .8093E-05 
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• 3473E- 04 .38166-04 
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The next page prints out integral properties, area and mass flux. It also prints out the 
iteration history of scalar potential equation at streamwise station 3. The ITER is the iteration 
number, PMAX is the maximum value of scalar potential, DELP is the maximum of change of scalar 
potential at one iteration. RHSMAX is the maximum residual of equation. The PTEST and RHS test 
are the measures of criterion of convergence test of equation. In the present calculation, when 
PTEST < 0.1, RHS < 0.1, and if RHMAX/RHMAXI < 0.005 (RHMAXI is the maximum residue at ITER=1) , 


the equation converges. 
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The next page prints out iteration history o£ stream function - vorticity equations at streamwise 
station 3 The ITER is the iteration number. RHO is the parameter for time scaling, MODEF and 
parameters relating to time scaling. DFF, PXF are the maximum value of charge of streamwise vorticity and 
maximum value of streamwise vorticity respectively. DPS, PXS are the maximum absolute value of charge of 
stream-function and maximum value of stream-function respectively. TS and TF are the measures of criterion 
of convergence test of equation. When TS < 0.1 and TF < 0.1 simultaneously ICONV value increases by 1. 

If IC0NV=5, then equation converges. 
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The following pages print out the Cartesian vertical coordinate ’Y', Cartesian spanwise 
coordinate *Z\ streamwise velocity ’u\ the vertical velocity ’v\ 
the spanwise velocity 'w* , the streamwise vorticity ’vor-x’ at grid point IY, IZ. 



3-rn s r a 1 1 o v at 

c?:nt-:ki. if u location c .0 , .0 


STATION 3 


* ** * 


Y 


12 = 1 2 3 4 

2 = • 1000K + 01 .20006+01 .30nG6+01 .4 O 0 06 + 0 1 


IY 

1 

Y 

•1000 K +01 

•5243C-0 1 

2 

•20nor +01 

•52486-D1 

3 

• 3000EZ + 01 

.52546-01 

4 

.<rOouF:+oi 

.52£24>D1 

5 

• S 00 Ci*. + 01 

.5272C-01 

6 

• 60001* +01 

♦ 5284 ”-0 1 

7 

• 7cuoi:+oi 

.52996-01 

8 

• 6000 £ + C 1 

•53186-01 

9 

• 9000 tl + 01 

.5342E-J31 

10 

• ioooe:+02 

•53735-01 

11 

• 1100t* + 02 

.54105-01 

12 

• ll'OOt +02 

• 5 4 5 8 r - 0 1 

13 

.13001+02 

.55175-01 

14 

.1 aoo-:*+o2 

.55516-01 

15 

.ltOO:>C2 

.56855-01 

16 

. if.ooe:+02 

♦5803F-01 

Kl7 

. 1700E + 02 

.59516-01 

18 

• 1 POOl+02 

.61*386-01 

19 

•1900L+0? 

♦ 6375K -0 1 

20 

.2000- -G2 

• £» 6 7 li * — P 1 

21 

.2100.*. + 02 

• 7067' -f l 

22 

.22003-02 

*75456-01 

23 

.230 OK +02 

*81726-01 

24 

.24003+02 

.8983' -01 

25 

• 2 tOOr. +02 

.1004-. + 00 

26 

.26001* +02 

•11456* GO 

27 

. 2700:. + 02 

.1333:: + CO 

28 

.21X03+02 

.16926+0 0 

29 

.2<r00;-. -02 

.l c 5r-6 + CG 

30 

.3roo r *02 

.250 0.:.- CO 


•52436-01 

•52486-01 

• 52545-01. 
.52626-01 
*52726-01 
.52845-01 
.5.2595-01 
•53185-01 
•53425-01 
.5 3735-01 
•54105-01 
.54575-01 
•55175-01 
•55515-01 

• 5 68 55-01 
•5 o 035-01 
.55515-01 
. f 1385-01 

• r- 3755-0 1 
•6 075 5-01 
. 7 0 5 f - 0 1 
.75 1 45-01 
•81715-01 

• 8^ci2c.-Ul 

• 1004/. >00 

. 1 1445+00 

• 1 333 5 + 0 0 
. 15526+^0 
•l r trl+OG 

• 2 4 : 5 + u 0 


•52435-01 

•52485-01 

• 5 2545-01 
•52626-01 
•52725-01 
•528*5-01 
.5299E-01 
•53185-01 

• 53425-01. 
•53735-01 
.54105-01 
.54575-01 
•55175-01 

• 5 5? 15-01 
.56855-01 
.58035-01 
.59515-01 
•6138E-01 
.63755-01 
•66755-01 
.70565-01 
.75445-01 
•41715-01 
•89825-01 
•1004C-00 
. 1 1*45+00 
♦13335+00 
. 15525+00 
.165 tOGC 
. 2 4 •> + r. - 0 0 


•52435-01 

•52485-01 

.52545-01 

.52625-01 

•52726-01 

.52645-01 

.52995-01 

.53185-01 

.53435-01 

.53736-01 

.54106-01 

•54585-01 

♦ 55 175 - C 1 
♦55925-01 
.56555-01 
. 58035-0 1 
.59516-01 
.61396-01 
.63756-01 
•667o r -01 

. 705*6-01 

. 75466-01 

♦ 81 / 3 - G 1 
.89856-01 
.10056+00 
♦11456+00 
.13336+00 
- 15925+00 

♦ 1 9 o 9 ;\ + G 0 
.2501 5*00 


.300006+01 

, .30000 6* Cl 


STEP SIZE •100005+C1 


•500C6+01 .60 6.0 L + 01 


.52435-01 .52436-01 

•52486-01 .52486-01 

.52546-01. . • 5254E-Q 1 
.52625-01. .5262E-01 

.52726-01 .5272^-01 

*52846-01 . 5284 6-01 

• 52=?9E-01 452996-01 

• 531 9F.-01 *53196-01 

.53436-01 *53436-01 

.53735-01 .53736-0 1 

• 54 1 IE- 01 .54115-01 

• 545 85-01 .5456 r. - Q1 

• 55175-01 • 55176-0 1 

.55926-01 .55926-01 

•56866-01 .56865-01 

*56046-01 . 58045-01 

.59525-01 .59535-01 

• 61405-01 .'614 05-01 

.63775-01 .63775-01 

• 66735-01 4 667 T 5-0 1 

• 70605-01. . 7061 r -01 

.754-5-01 .75466-01 

.81785-01 *31785-01 

.89915-01 • .89915-01 
.10056+00 .10055+00 

•11465+00 .11*65+00 

.13355+00 .13355+00 

.15^-6+00 •I5 r 45+0Q 

•19417 + 00 .19615+00 

• 250 * ~ + 0 u .25045*00 


7 8 

.70006+01 .80006+01 


.52435-01 .52425-01 

.52486-01 .52475-01 

.52545-01 .52535-01 

.52626-01 .52616-01 

*52725-01 .52705-01 

.52846-01 .52626-01 

• 5299 6r 0 1 .52976-01 

• 5318 E-,0 1 .53166-01 

.53426-01 .53406-01 

.53726-01 .53695-01 

.54106-01 .54066-01 

.54576-01 .54526-01 

.55166-01 .55106-01 

.55906-01 .55836-01 

.56846-01 .56756-01 

.58015-01 .57906-01 

.59496-01 .59356-01 

.61356-01 .61136-01 

.63716-01 .63496-01 

.^6706-01 .66436-01 

• 7 051 5-01 .70165-01 

.75375-01 .74936-01 

.01626-01 .81065-01 

.89705-01 .89005-01 

.10035+01 .99376-01 

•11*25+01 .11315+00 

.13305+01 *13155+00 

.15685+01 .15685+00 

.15545+01 .19275+00 

•24935+01 .24576+00 


q 

.90006+01 


.52346-01 

.52386-01 

.52446-01 

.52515-01' 

.52605-01 

.52726-01 

.52866-01 

.53046-01 

.53276-01 

.53555-01 

.53905-01 

.54346-01 

•5490E-01 

.55595-01 

.56476-01 

.57576-0 1 

.56965-01 

.60716-01 

.62325-01 

.65725-01 

.68295-01 

•73865-01 

.79725-01 

.87316-01 

.97225-01 

.11035+00 

.12805+00 

.15226+00 

♦16645+00 

• 237 It' ♦ C 0 


10 

.10006+02 


.51826-01 
.51866-01 
.51926-01 
'• 5198E-0 1 
.52076-01 
.52176-01 
.52316-01 
• 5247E-0 1 
•5268E-0 1 
.52946-01 
.53276-01 
♦5368E-0 1 
-.54196-01 
.54836-01 
.55646-01 
.56666-01 
.57356-01 
.59566-01 
.61616-01 
.64216-01 
.67516-01 
.71736-01 
.77156-01 
.•84176-01 
•9335E-0 1 
.10556+00 
• 121 8E + 0 0 
.14425+00 
.17596+00 
.22276+00 



Y Y 


1 

• 1 000E+01 

»50D7tl-D 1 

2 

•2000E+01 

• 5 0JL1E-D 1 • 

3 

• 3 OOOE + Ol 

.5016E-D1 

4 

•4000E+01 

•5022E-D 1 

5 

• 5 000 E + 0 1 

.5029E-D1 

6 

• 6 000 E + 0 1 

.5039E-DX 

7 

•7000E+01 

•5050E-C 1 

8 

• 8000E +0 1 

. 5063E-D 1 

9 

•9000E+01 

.soB^E-oi 

0 

• 1P00E + 02 

. 5 1 D ft “0 1 

1 

•1100 E+02 

.sisec-Gi 

2 

• 1 200 E+02 

.5173E-D1 

3 

•1300E+02 

.5219E-0 1 

4 

• 1 400 E+02 

• 5276E-D 1 

5 

•1500E+02 

*53^ 9E“0 1 

6 

• 1 COOE + 02 

•5439E-B1 

7 

•1700E+02 

.555^1 T.-Dl 

8 

• 1800E+02 

.5699F-D1 

9 . 

•1900E+02 

« 5Gfi2c. - 0 1 

0 

• 2 000 E+02 

*6 113F.-0 1 

1 

. *2 100 E + 02 

.64D8F-D1 

•2 

•2200E+0 2 

• 67GSE-0 1 

3 

•2300E+D2 

.7270F-D1 

4 

• 2400E+C2 

,7e97E-Dl 

5 

•2500E+02 

•S716E-01 

•6 

• 2600 E + 0 2 

.9G00E-01 

•7 

•2700E+02 

• 1 126E + 0 0 

8 

• 2.800E + 02 

.X326E»D0 

9 

•2900E+C2 

• 1 6 09E + D 0 

0 

. 3 OOOE+02 

• 2C27E *0 0 


12 = 21 
2 = • 21 0 OE + 02 


•468SE-01 .4052E-01 

•46B6E-01 • • 4 055E- 0 1 

• 4 69DE*'01 *4059E-D1 

• 4 695E-Q 1 .4063E-01 

•4702E-01 « 4 068E- 0 1 

• 4 71 0 E-0 1 .4075E-01 

•4720E-01 • 4 08 3E- 0 1 

•4732E-01 4 4 093E-01 

.4748E-01 ;4106E-01 

.4768E-01 .4123E-01 

• 4 793E-0 1 .4143E-01 

• 4 824 E^O 1 • 4 169E- 0 1 

• 4 864 E-0 1 « 4 201E-0 1 

• 4 913E-0 1 ^ 4 24 2 E-0 1 

• 4 97 5E-0 1 *4 292E-0 1 

• 5 0 5 3 E - 0 1 .4 356E-01 

•5151E-01 • 4 43 7 E- 0 1 

•527SE-0 1 • 4 539 E-0 1 

« 5 431E-0 1 .4667E-01 

• 5630 E-0 1 • 4 830E-0 1 

•5633E-01 .5037E-01 

•6205E-01 .5302E-01 

• 6620E-0 1 • 5643E-01 

• 7 157E-0 1 • 6 084 £-01 

• 7 65 9 E-0 1 • 6 659 E-01 

• P736E-0 1 .7421E-01 

.10035+00 • 8444E-C1 

• 1 175E + 00- *"9850E-0 1 

• 1417E+-00 • 1 1P4E + 0 0 

• 1 775E + 0 C *1476E+00 


22 23 

•2200E+G2 .2300E+02 


• 3214 E-01 
.3217E-01 
•3219E-01 
♦ 3222E-0 1 
•3226E-01 

* 323 2 E- 0 1 

• 32 38 E -01 
.32465-01 
• 3256E-0 1 

* 326 BE “0 1 
, 3284 E-01 

♦ 33 04 E-01 
♦332SE-01 
• 3359E-0 1 
• 3398E-01 

• 3447 E-01 
. 3509 E-01 
•35872-01 
. 3685E -0 1 

* 38 1 0 E-01 
•3965E-01 
•4172E-01 
•4432E-01 
•4770E-01 
•521 1 E-01 

* 5794 E-01 
•6577E-01 

• 7654 E- 0 1 
•9178E-01 
• 1 1 43E+00 


24 

.24005*02 


♦2222E-01 

• 2223E-01 
.2225E-01 

• 2227E-0 1 

• 2230E-01 
•2233E-01 
•2 237E-01 
•2243E-01 
•2250E-01 


• 2 25 BE - 0 1 
• 2269E- 01 

• 228 2E- 0 1 

• 2299 E- 0 1 
•2320E-01 

• 234 6E- 0 1 
• 2379E-01 
.24212-01 
•2474 E-01 

• 254 0 E-01 
•2625E-01 
•2732E-01 
•2870E-01 
•3047E-01 
•3275E-01 
•3574E-01 
•3969E-01 
.44992-01 


•5229E-01 
•626 1 E- 01 
•7787E-01 


25 

•2500E+02 


•1133E-01 

*1134E-01 

• 1 135E-0 1 
.1136E-01 
•1137E-01 

• 1 1 39E-0 1 
•1141E-01 
•1144E-01 

• 1 14 8E -0 1 

• 1 15 2E-0 1 
. 1157E-01 

• 1 164E-0 1 
4 1 172E-01 
•11B3~-01 
ill96?:-01 
41213^-01 
.1234E-01 

• 1 26 1 E-0 1 
•1294E-01 
•1337E-01 
41391E-01 

• 1461E-01 
• 1550E-0 1 

• 1.665E-01 
•1816~-01 
42016 E - 0 1 
*228 4 E-01 
•2652E-01 
43173E-01 
.3944^-01 


• 21 28 E-l 5 

• 2 1 26 E-l 5 
• 2123E-1 5 

• 2120 E-l 5 

• 2 1 1 6 E - 1 5 

• 21 10 E-l 5 

• 21 04 E-l 5 

• 20 96 E-l 5 

• 20 86 E-l 5 

• 20 73 E-l 5 

• 20 57 E-l 5 

• 2037 E-l 5 
•20 11 E-l 5 
•1979 E-l 5 

• 1940 E-l 5 

• 1 889 E-l 5 
• 1826E-1 5 
« 1746E-1 5 
« 1646 E-l 5 
•1518E-15 
♦ 1355E-1 5 

• 1 147 E-l 5 

• 8800 E-l o 

• 5342 E-l a 

• 8236E>17 
-.5152E-1 S 
-* 1318E-1 5 
-•2422 E- 15 
-•39 8 3 E-l 5 
-.6291E-15 


-•1133E-01 
- *1 134 E -01 
-•1135E -01 
-•1I36E-01 
-•1I37E-01 
-•1139E-01 
- • 1141 E -01 
-41144 E-01 
-•11 48 E-01 
-.1152E-01 

- • 1 1 57 E - 01 
-• 1.164 E-01 
-•1172E-01 
-• 1 185E -01 
-«rl 1 96E - 01 
-.1213E-01 
- • 1234 E -01 
-*1 261E -01 
- • 1 294 E- 01 

- • 1337 E-01 
-•1391 E -01 
-•1461E-01 
-• 1550 E -01 
-• 1665E -01 
-.1816E-01 
-•2016E-01 
-.22845-01 
-•2652E-01 
-.3173E-01 
-♦3944 E -01 


-• 2222E-0 1 
-• 2223E-0 1 
*• 2225E-0 1 
•• 2227E-0 1 
-•2230E-01 
-.2233E-01 

- • 2237E-0 1 
-.2243E-01 
• • 225 QE-0 1 

- .225 8E-0 1 

- • 2269E-0 1 
-.2282E-01 

- • 229 9E-0 1 
-♦2320E-01 

- • 234 6E-0 1 

- • 2379E-0 1 

- • 242 1E-0 1 
-•2474 E-01 
-•2540E-01 

- • 2 62 5E-0 1 
-• 2732E-0 1 

- • 267 OE-0 1 
-• 3047E-0 1 

- • 3275E-0 1 
-.3574E-01 
-• 3969E-0 1 

- • 44.99E-0 1 

- • 5229E-0 1- 
-•6261E-01 

- • 778-7E-0 1 


-.3214E-01 

- •3217E-0 1. 
-.3219E-01 
-•3222E-01 
-•3226E-01 
-•3232E-0 1 
-•3238E-0 1 
-.3246E-01 
- • 3256E-0 1 
- • 3268E-0 1 
- • 3284 E-0 1 
-•3304E-01 
-•3328E-0 1 
- •3359E-0 1 
- •3398E-0 1 
-•3447E-0 1 
- • 3509E-0 1 
-• 3587E-0 1 
-•3685E-01 
-•3810E-01 
- • 3969E-0 1 
-•4172E-0 1 
- • 4432E-0 1 

- • 477 0£-0 1 
- .5211E-0 1 
-•5794 E-01 
-•6577E-D 1 
-•7654E-01 
-•9178E-01 

- • 1 14 3E+0 0 


2 6 27 28 2 9 3 0 

• 26 0.0 F+ 0 2 .2700E+D2 .2800E+02 .2900E+02 *3000E+02 



Y 

• 1 000 E+Cl - • 4 0 52E-0 1 -.4683E-01 

• 2 OOOE + 0 1 -.4055E-D1 .-.4666E-01 
.3 OOOE+0 1 - • 4 059E-C 1 -.4690E-01 
.40005+01 -.4C63F-01 -•4C95E-01 

• 5CGCE + 0 1 -.4066E-DI -.47025-01 

• 6 OGO E + 0 1 -.4075 f. - Cl -•*710001 

• 7000 E + 01 - * 4 0 6‘3 E .- 0 1 -.*720001 
*60005*01 — *40936-01 -.47325-01 
.90005+01 - * 4 106E-0 1 -.47485-01 
•1000E+02 - • 4 1235-C 1 -.47685-01 
*11005+02 - . 4 1 43E-0 1 -.4793E-01 
• 12005+02 - • 4 1 69E-C 1 -.4E24E-01 
•1200 E+02 -.42016-01 -* '-6645-01 

• 1 400E + 02 -.42426-01 -.* n 13E-ni 

• 1500E + 02 -.42921. - 01 -.49755-01 

• 1 600E+02 -.4356 E-01 -.50535-01 

• 1 700E + 02 - • 4 437E-D 1 -.51515-01 
•1800E+02 -.4539 E-01 -.5275E-01 
• 1900E+02 -.4667E-01 -.5431E-01 

• 2 000 E + 02 -.40305-01 -.5630E-C1 
*21005+02 -.5037;-Cl -.58B3E-01 


.22005+02 

•23005+G2 

.24005+02 

•2500E+02 

• 2 CO 0 E + 0 2 
• 2700E+02 
.2TO0E+O2 

• 2 c 00E + 02 

• 3 COPE +02 


• 5 30 2 r - 0.1 

• 5643E -0 1 
• 6084E- 0 1 
•6659E-D1 
♦7421E-01 
•8444E-D1 
•9S50E-D1 
.1104E+OO 
•1478 E *00 


-.62056-01 
-• 6620 E-01 
• • 7 157E-0 1 
-.7859E-01 

- • 8786E-0 1 

- • 1 003E + 00 
“• 1 175E + 0 0 
-.1417E+00 
-•1775E+00 


-;5007E-01 -.5182E-01 
-•5011E-01 -*51865-01 
-• 5016E-0 1 - • 51 92E-01 
-.5022E-01 - • 5198F - 01 
- * 5 02 9 E- 0 1 -.5207F-01 
-.5039E-01 “.5217E-01 
-.50 5 OF- 01 523 1 E-01 

- « 5 065E- 0 1 - * 52 4 7 E-0 1 
- i 5 084E- 0 1 - • 5268E- 01 
-♦5107E-01 - • 52 94 E- 01 
- . 5 136E-0 1 -.5327 E-01 
-.5 17 3 E-01 -.536 8 E-01 
-.521^ E-01 -.54 19E-01 
- . 5 276E- 0 1 -.5483E-01 
-•5349E-01 -.5564E-01 
-•5439 E-01 -♦ 56 66 E-0 1 
- . 5 554 E-0 1 - • 5795E- 0 1 
- . 5 699E-0 1 -.5956E-01 
-.5882E-01 -.61615-01 
-•6113E-01 - • 64 2 IF- 0 1 
- • 6 408E-01 -.6751E-01 
-• 67F5E-01 -.71736-01 
- . 7 270E-0 1 -.7715E-01 
-.7O97E-01 -• 64 1 7 E-01 
- • 8 716E-0 1 - * 9335E-0 1 
-•9800E-01 -.1055E+00 
- . 1 126E+ 0 0 - . 12 1 6E + 00 
-•1326E + 0C - • 14 4 ?E + 00 
1609E+00 -.175^5+00 
-.2027E+OC -♦2227E+00 


•5234E-01 -.5242E-01 

• 5238E-01 - ; 5 24 7E-Q 1 

• 524 4 E- 01 -♦ 5253E-0 1 
♦5251E-01 - • 5261E-0 1 
•5260E-01 -.5270^-01 

• 5272E-01 - 4 5 2 8 2 “ - 0 1 

.5266E-03 -.5297=3-01 

.5304 E-01 5316E-01 

•5327E-01 -.5340E-01 
•5335E-01 - • 5389E-0 1 

• 5390E-01 - « 5406E-0 1 

• 54 34 E-01 -*545? r .-01 
•5490E-01 - • 551 OE- 0 1 

• 555 9E- 0 1 - 4 55P 3E- 0 1 
♦5647E-01 - ♦ 5675E-0 1 

• 5757 E- 0 1 -.5790E-01 
• 5896E- 01 -.593 5 E-01 
•6071E- 01 -.6118E-01 

• 6292E-01 -*634 q E-0 1 
•6572E-01 -.6643E-01 

• 6529E-01 - . 7 016E-0 1 

• 7366F-01 - « 7493E-0 1 
.797 2E- 01 -.8106 r -01 

• 873 IE- 0 1 -♦890 0 E-01 
• 9722E-0 1 -.9937E-01 

• 1 1 03E + 00 -;1131E + 00 

• 1230E+ 00 -♦ 1315E + 00 
•1522E+00 -• 156FE+00 
•1K64E+00 -.1927E+00 
•2371E+00 -• 2457E+00 


-• 5243 E-0 1 -•5243E-01 
- • 5 2 4 8 E - 0 1 -•5248E-01 
- *5254 E-0 l -.5254E-01 
-.526? E-0 l -.52 62 E-01 
- • 5272 E-0 1 -.5272E-01 
-•5284 E -01 -.5284E-01 
- . 5299 E-01 - .5299 E - Cl 
-.5318E-01 -.5319E-01 
-.5342 E-01 -.5343E-01 
-•5372E-01 -.5373E-C1 
-•5410 E-01 -.5411E-01 
-.5457 E-01 - • 5458 E -01 
5516E-0 1 -•5517E-01 
-•5590 E-0 l - • 5592 E-01 
-.5684 E-01 -.5686E-01 
-• 58 0 1 E-0 l - *58 04 E -01 
- » 5949 E-0 1 -.5953E-01 
-•6135 £-0 l - *61 40 E -01 
-•6371 E-0 l -.6377E-01 
-.6670 E-01 - •6678E - 01 
-.7051 E-0 L - .70 61 E -01 
-.7537 E-0 l -.7549E-01 
-.8 162 E- 01 - • 81 78 E - 01 
-•8970 E-01 -.8991E-01 
-.1003 E+ 03 -.1005E+00 
-• 1 1 42 E + 0 3 - • 1 1 46E +00 
-•1330E+03 -*13355+00 
-•1588E+03 -.1594E+00 
1*954 E + 03 -*19615 + 00 
- . 24 93 E + 0 3 - • 25 04 E + 00 


• 5243E-0 1 -.5243E-01 
•5248E-01 - • 5248E-0 1 

• 525*4E“01 -.5254E-01 

• 5262E-0 1 - • 5262E-0 1 
•5272E-01 -•5272E-01 

• 5284 E-0 1 -.5284E-01 
. 5299E-C 1 - *5299E-0 1 
• 5319E-0 1 - • 5318E-0 1 

• 5 34 3E-0 1 - • 534 3E-0 1 
.5373E-01 -•5373E-01 
•5411E-01 -•5410E-01 
.5458E-01 -.5458E-01 
.5517E-01 -.5517E-01 
•5592E-01 -.5592E-01 

• 568 6E-0 1 - • 5685E-0 1 
. 580 4 E-0 1 - *5803E-0 1 
• 5952E-0 1 -•5951 E-01 
.6140E-01 -• 6139E-0 1 

• 6377E-0 1 - * 6375E-0 1 
•6678E-01 -.6676E-01 

• 7 0 6 OE - 0 1 -WOSSE-Ol 
•7549E-01 - • 7546E-0 1 
•8178E-01 -.8173E-01 
•B991E-01 -.8985E-01 
•1005E+00 -.1005E+00 

• 1 146E + 0 0 -•1145E + 00 

• 1 335E + 0 0 -• 1333E + 0 0 

• 159 4E + 00 - * 1592E + 0 0 
•1961E+00 -.1959E+00 

• 250 4£^0 0 -*2501E + 00 



• 1 000 E + 01 
•2000E+01 

• 3 00GE + 0 1 

• 4 000E + 01 
•5oooe+oi 
•6000E+01 

• 7 000E + 01 
•8000E+04 
•9Q00E+0L 

• 1 000E + 02S 

• IIOOE+OSK 

• 1 200E+02 

• 1 300 E + 0 2 
» 1 400E + 02 

• 1 500E +02 

• 1 600E + 02 

• 1 70QE + 02 
* 1800E+02 

• 1900 E + 02 
• 2000E+02 
• 2100E+02 
• 2200E+02 
•2300E+02 

• 2 400 E + 02 
•2500E+02 
•2600E+02 

• 2700E +02 

• 2800E + 02 
►2900E+02 
■ 3 OOOE + 02 


“ • 5243E-0 1 
-♦5248E-0 1 

- • 5254E.-0 1 
-.5262F-01 
-.5272E-G1 
-•5284E-01 
-•5299E-0 1 
-.5318E-01 
" *5 342E -JO 1 

■>- • §*3 7 3 E - 0 1 
35* 541 OE -01 
5457 E -01 
-.551 7E-0 1 
-• 559 1E-0 1 
-.5685E-01 
- • 5803E -0 1 
-.59S1E-01 
- • 6 1 38E -0 1 
- « 6375E -01 
- • 6675E-D 1 
-•7056 E -01 
-•7544E-01 
-*817 IE-01 

- • 8982E-0 1 
-.1004E+00 
-.1144E+00 
-•1333 E + 00 
- • 1 592E + 0 0 

- • 1958E +0 0 
-•2499E+00 


-•5243E-01 

- • - 248E-0 1 
-.5254E-01 
-•5 26 2 E -01 
“ • 5272E-0 1 

• - • 5284E-01 
-• 5299E-0 1 
- • 5318E-0 1 
-•5342E-01 
-•5373E-01 
-• 5410E-0 1 
“ *54575-0 2 
-•5517E-01 
-.5591E-01 

- • 5685E-0 1 
-.5803E-01 

- • 5 95 IE- 0 1 

- • 6 138E-0 1 
* • 6 375E-0 1 

- • 8 675E-0 1 
-• 7056E-0 1 
-•7544E-01 
-• 8 171E-0 1 

- • 8 982E-0 1 
-•100.4E + 00 
- • 1 14 4E + 00 
- • 1 333 E + 0 0 
-•1 592E+00 
-• 1 958E+00 
-•2499E+00 


-•5243E-01 
-.5248E-01 
- • 5254E-01 
-.5262E-01 
-•5272E-01 
-.5284E-01 
-•5299E-01 
-.5318E-01 

- ♦ S 342E- 0 1 

- ^ 5 373E-0 1 

- • 541 Ot-Ol 
-.5458E-01 
-.5517E-01 
- • 5591E-01 
- • 5685E-01 
-•5803E-01 
- • 5 951 E- 0 1 

- • 6 138 E -0 1 
- • 6 375E-0 1 
- • 6 675E-01 
-•7057E-01 

- • 7 545E-01 
-•8172E-01 
- • 8 983E-01 
-•1OO4E+0Q 
- • 1 145E + 00 
-*1333c+00 
- • 1592E + 0 0 
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•2400E+02 -.1730E+00 
•2500E+02 -• 1728E + D 0 
. 2600E + 02 -..1725E+00 
.2700E+02 -.1722E+00 ■ 
•2800E+02 -.1717E+0Q 
•2S00E+02 -.1710E+00 ■ 
.3000E+02 -.1700E+00 ■ 


• 32 0 0 E <• 0 
- • 2 0“6E + 0 0 
- • 2 096E*0 0 
-.20965+00 
-.20965+00 
- .2 096E + 0 0 
-.2096E+00 
- . 2 096E + 0 0 
-.20965+00 
-.2096E+0Q 
-.2096E+00 
-.2096E+00 
-.2096E+00 
-.2096E+00 
-.2096E+00 
-.2096E+00 
- • 2 096E + 0 C 
- . 2 096E + 0 0 
- . 2 096E + 0 0 
- . 2 096E + 0 0 
-.2096E+00 
- . 2096E + 0 0 
-.2096E+00 
- . 2 096E + 0 0 
-.2097E+00 
-.2097E+00 
-.2097E+00 
- . 2 097E + 0 0 
.2098E+00 
. 2 0V9E + 0 0 
• 2100E + 00 


2 . 3300E+02 

-•2500E+00 
-.2500E+00 
- *2 SO 0 L + 0 0 
-.25005+00 
-.2500E+00 
-•2500E+00 
-.2500E+00 
-•2500E+00 
-.2500E+00 
-.2500E+00 
-.2500E+00 
- .2500E+00 
-.2500E+00 
-.2500E+00 
-.2500E+00 
-.'2500E + 00 
-.2500E+00 
-.2500E+00 
-.2500E+00 
- • 2500E+00 
-.2500E+00 
-.2500E+00 
-.2500E+00 
-.2500E+00 
-.2500E+00 
-.2500E+00 
-.2500E+00 
-i2500E+00 
-.25005+00 
-.25005+00 
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T7” 

1 

0 

3 

4 


7 = 

.10001+01 

.20006 + 01 

.3ooo:->oi 

.40005+01 

rr 

Y 





l 

.10006+01 

.78485-1 1 

.0 

.0 

• 0 

2 

.20005+01 

.585+1 -01 

. 5 <: 6 2 K ”01 

.bti72E-01 

. 59 ft 65 - 0 1 

3 

.30005+01 

. 1 £?3>>G0 

. lC23c+00 

.1 025E-00 

*10515+00 

4 

.40005+01 

-14ia**-00 

♦ i ‘iloil+oo 

. 1 419E+00 

.14516+00 

5 

.50005+01 

, i*:ois>rc 

• 1 i * 0 3 5 + 0 ft 

. 1 603E+00 

« 1 S 3 7 r +00 

6 

.60005 + 01 

.220::*:+ 00 

.2: ?r>li.+OC 

• 2i;01E *00 

.2236:5+00 

7 

.70005+01 

+00 

.Ti^7L + C0 

.^£27E»00 

.26586+00 

8 

.80005+01 

.3ot3::+co 

.3 082 ’- + 00 

.30r(lE + CO 

•310KE+0C 

9 

.90005+01 

.35665+00 

. 3 565L. + 0 C 

• 3 563E + 0 0 

.35886+00 

10 

. 1000 + 02 

. 4 06{;L' + 0 0 

+4CC 7C+01 

• 4 066E + C0 

.40886*00 

11 

. 13005+02 

• ft 5 7 9 -* + 0 0 

r457cU>0li 

• A 570E+ 0 0 

.45956+00 

12 

.12006+02 

#50835-60 

#508^+00 

• 5 OfJOE + DO 

. 50S6K+00 

13 

.13005+02 

.55645 + 0 0 

.566^6+00 

•556RE+00 

• 55 7 4 5 + 00 

14 

.14006+02 

. 6 0 1 1 r . + 0 0 

.60105+ OC 

.&00ft'i»00 

.60175+00 

15 

.15006+0? 

♦ 6413 r +00 

.64 126+00 

. 641 OE + OO 

♦6417-+00 

16 

.16005+02 

.676 8 V +00 

.C 7f,7r +0G 

. 6 766c +0 0 

. 677 Or +CC 

17 

.17006+02 

.70805*00 

•7uftOr +00 

♦ 7 07fcE*00 

.70 ft 15+00 

1 3 

♦18005+02 

•73585+00 

#73585+00 

. 7 3 c jf,L*CG 

*735? ; * + 0C 

13 

•190Dr+G2 

.76105+00 

• 7 £ 1 G r. + 0 C 

. 7 EOF E +00 

.76095+00 

20 

*20005+0 c 

♦ 7c; 4 61.4 0 u 

4 7 ; *\ f: + 0 ( 

. 7 5:45000 

• 7 8 4 5 ir +C0 

21 

. 2 1006 + 0? 

.P0‘>5i +00 

• f U\ 45+ Of 

• 1 0V3E ♦ OC 

. 80 9.31.* + 00 

22 

.22005+02 

-c 37 6{ +C0 

. ; 5 7 6 £ + (1 C 

.6375E+0& 

.83756+00 

23 

• 2700i>02 

. ft 6V6i. * fj u 

• 8 i i; 6 5 4 5 r 

♦ fc.t G EE + Oi. 

• s6‘ir..T + l , C 

24 

• 2 4 0 0 r +C2 

♦9061r *♦ C £ 

• c r-bic>no 

. vOECi +00 

.9061 r+CC 

25 

• 250c; : +02 

. f 4C6.. + C C 

+00 

• ^ r; 5 + II c 

• 9467*' +00 

26 

*2roo:+C:2 

• c ?. .? 4 * +C0 

# r / 2 ft 5 ♦ 0 (. 

.5tr6E+0C 

. hi'* 3 C 5 + 0 0 

27 

.2700- +02 

• 1131 a ‘ + 01 

. 1 1 : +^>1 

. 3 ID L + C 1 

• 1 1 1 c ? + 0 1 

2ft 

.2*006+02 

- UU<r -c 1 

• 1 115.- +01 

. 1 1 r- 7 : . + 0 1 

.1057" +01 

23 

• 2 c G U r + 0 ? 

.10* 1V.-01 

. 3 08 7. : 4 Cl 

. 1 0.-2E + U1 

• 1 0 7 3 > + fi 1 

30 

. 3 C 0 C • +05 

. 1 flt 5; ! -rf 1 

. 3 GbGr. + 01 

. 1 05 ft : + C 2 

. 104r "+G1 


5 

.50005+01 


.0 

.5*235-01 
•10415+00 
. 14446+00 
.18331+00 
.22356+GG 
- 26595+00 
•31126+00 
.35935+00 
.40955+00 
.46036+00 
*51046+00 
.35825+00 
•60236+00 
•64211+00 
•67735+00 
•7083 6 +G0 

• 7358*:: + 00 
.76096+00 
.7845* +00 
•80936+00 
.83765+00 
. r697f+00 
•90635+00 

• c * 43 5 5 + 00 
# 96366+00 
• 11005+01 
.1087-* + 01 
. 10645* 01 
.103 c t + C 1 


6 

•6000F+01 


7 e 

.700 06+01 .80006 + 01 


9 10 

.90006+01 .10006+02 


;o 

• 584 25-0 1 
. 1 033 r + 0 0 
.14329+00 

• 1 P 1 7T + 00 
.22149+00 
.2635^+00 
k3085i+00 
•35666*00 
. 4067-+00 

• 45776 + 0 0 
.50809+00 
.5 C 62~+C0 
♦6007^+00 
.6408^+00 
. 67639+ CO 
.70769+00 
.7357- +00 
.76059+00 

• 784 2 r + 0 C 
. 8P919-0C 
. 8374 "+0C 

• B6 r '5" + G0 
•9C6?~+nC 

• 946* *>0C 
. 984 1 '’♦CO 
. 1 0879 + 0 t 
.10759+01 
.1054 >01 

• io? f *' + ni 


. o 

♦57796-01 
. 10299+00 
.1421 6+00 
.18019+00 
.21929+00 
.2610:! + 00 
.30579+00 
.35316+00 
.40309+00 
.45406+00 
.50459+00 
.55299+00 
•59786+00 
.83846+00 
.67449+00 
.70606+00 
.73409+00 
.7594 6+C 0 
.78316+00 
.8080 5 + 00 
■ 6 363 5 + 0 0 
•*6846+03 
.30515+00 
*94569+00 
•96405+00 
.10746+01 
#10635+01 
.10436+01 
.10215+01 


• 0 

.56816-01 

.10165+00 

• 14 136 +C0 
.17916+00 
.21776+00 

• 2567 E + 00 
.30265+00 
•34955+00 
.39896+00 
.44956+00 
.49996+00 
.54845+00 
.59376+00 
.63475+00 
•67116*00 ■ 
.70305+00 
.73146+00 
.75715*00 
•78085+00 
.80555+00 
.83356+00 
•86545+00 
.30195+00 

• °4 26 6 ♦ CO 
.9821 r + CC 
.10615+01 
.10505 +C1 
.10326+01 
.10155+01 


• 0 

.57506-01 
.10055+00 
.13365+00 
.17615+00 
.21416+00 
♦25445+00 
« 297 65 + 0 0 
.34425+00 
.39296+00 
.44316+00 
.49345+00 
.54206+00 
.58765+00 
.62905+00 
.66576+00 
.69816+00 
•726B5+00 
•75275+00 
.77655+00 

• S 00 95 + 00 

• t 26 46 + 0 0 
.859 85+00 
.89575+00 
.93615+00 
.97755+00 
•10465+01 
.10385+01 
.10235+01 
.10105+01 


.58746-01 
• 996 1 E-0 1 
.13616+00 
•17226+00 
.20976+00 
.24966+00 
•2925E+00 
•3386E+00 
.38736+00 
.43756+00 
• 4877E+Q 0 
.53636+00 
•58196+00 
.62346+00 
.66046+00 
.69306+00 
.72186+00 
.74796+00 
.77196+00 
.79606+00 
•8227E+00 
. 85346+0 0 
.88856+00 
.92846+00 
.97116+00 
.10386+01 
.10296+01 
.10186+01 
♦1006E+01 
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IZ* 

z = 


• 1 000 E + 0 1 
•2000E+01 

• 3 00 QE ♦ 0 1 

• 4 000 E + 0 1 
•5000E+01 
•6000E+01 
♦7000E+01 

• 8 00 0E + 01 
•9000E+01 
•1000E+02 
•1100E+02 
•1200E+02 

• 1 30 0E + 02 
•1400E+02 
* 1 50 0 E + 02 

• 1 600 E + 02 
• 1700E+02 

• 1 P00E ♦ 02 

• 1 900 E +0 2 

• 2000E + 02 

• 2 10GE + 02 
•2200E+82 
•2300E+02 
• 2400E+02 
*2 500 E + 02 
•2600E+02 

• 2700E + 02 
•2800E+G2 

• 2500 E + 02 

• 3 OOOE + 02 


IZ = 

z = 


•1000E+01 
♦2000E+0 1 

• 3 000E + 01 

• 4 000E + 01 

• 5 000 E + 0 1 
•6000E+01 

• 7 000E + 0 1 i 

• 8 000 E + 0'1 

• 9 0 0 0 E + 0 1 

• 1 000E+02 
• 1100E+02 

• 1 200.E + 0 2 

• 1300E + 0 2 
• 1 400E+02 
•1500E+02 
« 160QE+02 

• 1 70QE + 02 
•1800E+02 
• 1900E+02 
•2000E+02 

• 2 100E + 02 

• 2200 E + 02 
•2300E+02 

• 2 40 0E + 02 

• 2 500 E + 02 
•2600E+02 
•2700E+02 
•2800E+02 

• 2 500E + 02 

• 3 000E + 02 


II 


12 


13 


14 


•1100E+02 *12006+0? • 1300E+02 *1400E+02 .1500E+02 .16DL0E+02 .1700E+02 .1800E+02 •1900E+02 .2000E+02 


T6" 


T r 


"IS" 


19 


.0 

• 54 05E-J3 1 
•9138E-0 1 

• 1254E + 0 0 

• 160 9E + 0 0 

• 1 983E + D 0 
•2366E+00 
•2823E+00 

• 3293E + 0 0 
•3750E+00 
•4302E+00 
•4813E+D0 

• 5307E + 0 0 
•5771C+C0 
•6193 K +00 

• 656b6 +0 0 
.68996+00 

• 7 1 5 2 6 + E 0 
.74566+00 
.77006. + C0 

• 7'M3i>C 0 
.H2096+C0 
.35.13c ♦DO 
. 606 16 + D G 

• 9255c ♦0 0 
.960O6+CO 

• 1 029 £ ♦0 1 

• 1 0226 + C 1 
•1013E+8 1 

• 1 004E + D 1 


21 

.21001 + 02 


.44631-01 
. 7 6966-0 1 
.1073E>00 
.13846+00 
.170-5c + C0 
•20686+00 
•24706+00 
.29236+00 
.34291+0 0 
•39796+0 0 
. 45546 + 00 
.51306+00 
.56626+00 
• 61U66 + 0 0 
.66266+00 
•7005E+D0 
.73246+00 
.76006+00 
.78466+00 
.80936+00 
.63676+00 
.66776+00 
.90276+00 
i9413E+00 
.97946+00 
.10036+01 
.99936+00 
.99706+00 
.59366 + 0 0 


.0 

.66406-01 
.10796*00 
.14466+00 
.16416+00 
.22446+00 
.26736+00 
.31356+00 
.36216+00 
.4 1246+OC 
.46316+00 
.51246+00 
.55696+00 
.60166+00 
.63976+00 
.67336+00 
.70306+00 
.72956+00 
.75396+00 
.77696+00 
.60056+00 
.02706+00 
.05756+00 
.65226+00 
.63146*00 
• 572&C. + 00 
.10236+01 
.10166+01 
.10116+01 
.10026+01 


22 

.22006+02 


.0 

.30536-01 
.72396-01 
.10336+00 
•134B6+CG 
. 1 6605+00 
.204 6. 6 + 00 
.24566+00 
.29176+00 
.34296+00 
.35026+00 
.45566+00 
.51306+00 
.56706+00 
.61796+00 
• 6 61 6.6 + 0 0 
.65926+00 
.73096+00 
.75026+00 
.70276+00 
.00706+00 
.63376+00 
.06396+00 
•B90OE+OO 
.93616+00 
.97476+00 
.10046+01 
.997.16 + 00 
.995.46 + 00 
.55256+00 


.64206-01 
.11556+00 
.15866+00 
.20066+00 
.24406+00 
.26876+00 
.32556+00 
.38406+00 
.43346+00 
.40246+00 
.52566+00 
.57386+00 
.61406+00 
.65006+00 
.60106*00 
..71016 + 00 
.73566+00 
.75566+00 
.70246+00 

• f 0£ 2c + 0 0 
.33346+00 
.66466+00 
.90016+00 
.52966+00 

• 57 ‘3 56 + 00 
. 1 0176+01 
. 1 0146 + 01 
.10096+01 
.55906+00 


23 

. 2300E + 02 


.45306-01 
.30266-01 
. 11126 + 00 
.14246+00 
. 1 7596+00 
. 21256 + 00 
•25336+00 
. 2SeB6+00 
•34R9E+00 
.40256+00 
.45806+00 
.51346+00 
.56646+00 
^61506+00 
.65816+00 
.65536+00 
• 7 272E + 0 0 
•75506+00 
.70006+00 
•80456+00 
.83146*00 
. b6196+00 
.85636+00 
.‘93466 + 00 
.97366+00 
. 1 0036+01 
.99576+00 
.99336+00 
;99116+00 


.54356-01 
.11456+00 
.16246+00 
.20396+00 
.24.62^ + 00 
.28986+00 
.33556+00 
.38316+00 
.43236+00 
.48146+00 
. 52946 + 00 
.57476+00 
.61636+00 
.65346+00 
.68605+00 
.71496+00 
. 7 4 0 8 C ♦ 0 0 
.76446+00 
.78776+00 
.81166+00 
.83976+00 
.87156+00 
.90766+00 
.94725+00 
.98546+00 
.10155+01 
.10126+01 
.10076+01 
.99846+00 


24 

>2 4 0 OE 


►02 


.53016-01 
.91676-01 
.12646+00 
.16086+00 
» 15-626 *00 
.23466+00 
.27686+00 
.32306+00 
• .37316*00 
.42606+00 
.47986+00 
.53256+00 
.58236+00 
.62766+00 
.66766+00 
. 70 22E + 00 
.73226+00 
.75896+00 
.78316+00 
•80766+00 
.33516+00 
. 86626 + OD 
.90146+00 
•9402E+00 
.97836+00 
.10026+01 

• 9938E + 0 0 

• 99 HE* 00 
.9894E+0D 


. 0 

.65116-01 
.11386+00 
.16586+00 
.20806+00 
.25056+00 
.29376+00 
.33926+00 
.38696+00 
.4358C+00 
.48516+00 
.533 1E + 00 
.57846+00 
.61995+00 
.65696+00 
• 6 0 9 5 1 ♦ 0 0 
.71836*00 
.74416+00 
►76816+00 
.79106+00 
.81546+00 
.64366+00 
.87506+00 
.51216+00 
.55166+00 
.56966+00 
.10146+01 
.10105+01 
.10056+01 
.99786+00 


25 

.25006+02 


• 0 

.54376-01 
.96186-01 
.13366+00 
.17076+00 
.20906+00 
• .2-5006+00 
.29446+00 
.34216+00 
.39246+00 
.44456+00 
.49666+00 
.54696+00 
.59396+00 
.63635+00 
.67386+00 
.70656+00 
.73536+00 
.76136+00 
•7854 E+ 00 
.81046+00 
.83906+00 
.87135*00 
.90776+00 
•9473E+00 
•5841E+00 
.99366+00 
•9912E+00 
•9802E+00 
.98826+00 


.65325-01 
.11955+00 
.161 75+00 
42036-+00 
.24425+00 
•28645+00 
.33045+00 
.37735+00 
442665+00 
.47676+00 
.52616+00 
.57326+00 
.61665+00 
.65525+00 
.68925+00 

wigo^ + oo 

.74556+00 
. 76985*00 
.79295+00 
.81755+00 
.84605+00 
.07835+00 
.91485*00 
» 95425 + 0 0 
.99046+00 
.10125+01 
.10085+01 
.10035+01 
.99715+00 


26 ' 

• 26 OD 6+ 02 


.56935-01 

.10325+00 

.14205+00 

417935+00 

.217 7" +00 

. 258 85 + 00 

.30326+00 

.35075+00 

.40095+00 

445245+00 

.50365+00 

.55295+00 

.59885+00 

.64005+00 

.67655+00. 

470045+00 

.73675+00 

476245+00 

.78655+00 

481185+00 

•8410 r +00 

.87405+00 

.91135+00 

495165+00 

498755+00 

.99435+00 

.98965+00 

.98615+00 

.98665+00 


.0 

.64606-01 
.11456+00 
.15666+00 
.19666+00 
.23646+00 
.27776+00 
•3217E+00 
.36886+00 
.41836+00 
. 46S 1 6+00 
.51966+00 
.56816+00 
.61295+00 
.65305+00 
.68815+00 
.71866+00 
.74576+00 
.77036+00 
.79366+00 
.81846+00 
.84706+00 
.87946+00 
.91596+00 
.95525+00 
. 99096 + 00 
.10105+01 
•1007E+01 
. 1002E + 0 1 
.99625+00 


27 

.27006+02 


.0 

• 5834 6-0 1 
.10356+00 
•1433E+00 
.18116+00 
.2201 E+0 0 
.26186+00 
.30666+00 
.35426+00 
.40446+00 
.45586+00 
.50686+00 
. 5556 E+0 0 
.60086+00 
.64166+00 
.67766+00 
.70916+00 
.73726+00 
.76276+00 
•7867E+00 
.81236+00 
.84156+00 
.87486+00 
.91236+00 
•9529E+00 
.98876+00 
.98956+00 
.98636+00 
.9&44E+00 
.98506+00 


.59636-01 
.1039E+00 
.14716+00 
.18496+00 
•2239E+00 
*26406+00 
.30716+00 
.35396+00 
•4039E+00 
.45586+00 
.5079E+00 
.5584 E + 00 

• 6054E + 00 
•6475E +00 
.68436 +00 

• 7162E +00 

• 7442E +00 
.76946+00 
.79296+00 
.81786+00 
.84646+00 
.87885+00 
.91526+00 
.95446+00 
.99016+00 
.10096+01 
.10056+01 
.1002E+01 

• 9953E +00 


.62085-01 
.10186+00 
.13755+00 
.17235+00 
.20875+00 
.24706+00 
.28925+00 
.33566+00 
.36596+00 
.43876+00 
.49266+00 
.54545+00 
.59506+00 
.63976+00 
.67875+00 
.71236+00 
.74146+00 
.76735+00 
.75126+00 
. 81616+00 
.84466+00 
.87676+00 
.51295+00 
.952 16+00 
.98826+00 
.10876+01 
.1 O04E+O1 
.10816+01 
.99466+00 


28 29 

.28006+02 .29006+02 


.59396-01 

•1043E+00 

.14396+00 

. 18216+00 
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- « 1 1 8 b .: - 0 5 
.56225-0 2 
.11885-01 
.15335-01 
.24545.-01 
.30535-01 
.35365-0 1 
.395SE-01 
.4279E-0 1 
.45135-01 
.46645-01 
.47505-01 
.48215-01 
.48695-01 
.40-295-0' 1 
.50535-01 
.53456-01 
.59185-01 
.68205-0 1 
.74805-01 
.74425-01 
. 7 057: - 0 1 
.63305-01 


.21685-03 

- .72645-02 
-.1170 -- 01 
-. 1 4365-01 
-.1(113-01 
-.17151.-01 
-.17715-01 
-. 1 1 0 35-U 1 
-.16315-01 
-.16715-01 
-.1^526-01 

- . 2 C G 3 1*. - 0 1 
-. 20755-01 
- • 2 12o6-0 1 
-.21346-01 
-.20685-01 
-. 1522E-01 

] 6895-01 
-.1 3.765-01 
- . c 6635-02 
+ .53845-02 
-.60253-03 
.47416-02 
.11165-01 
.15575-01 
.31146-01 
.*3056-01 
.50135-01- 
.54036-01 
. 54426-01 


, 1 3945-03 
.75735-02 
.12365-01 
15765-01 
• 1 8675-01 
.21 '-06 - 01 
.25185-01 
.28825-01 
.32985-01 
.37696-01 
4 2855-01 
.4 8246-01 
.53475-01 
.56215-01 
.61845-01 
.64246-01 
.64935-01 
.63935-01 
.61175-01 
.56866-01 
.51215-01 
.44506-01 
-.36646-01 
-.27975-01 
-. 1 749E-01 
4 34 1 6-02 
.11266-01 
. 2.4815-01 
.36806-01 
. 43096-01 


.48875-04 
-.88566-02 
■ . 1455”-01 
■ . 19025-01 
.23295-01 
.27 7 3r -01 
.3261 5 -Cl 
.38166-01 
.44456-01 
.51 476-01 
.59055-01 
668* -.-01 
74386-01 
.81225-01 
.86906-01 
.90915-01 
-.93165-01 
-.93515-01 
-.. 91 8 C 5 - 0 1 
-.88526-01 
-.83596-01 
-.77285-01 
-.69475-01 
-.60165-01 
-.46985-01 
-.35375-01 
-. 18365-01 
— .6016'; — 03 
.17426-01 
.31436-01 


.10396-04 

-.90086-02 

- .15006-01 
-.20146-01 

- .25205-01 
-.30585-01 
-.36.536-01 
-.43225-01 

- .50686-01 
-.56346-01 
-.6747- -01 
-.76225-01 
-.64715-01 
-. c 2355 -01 
-.98935-01 
-.10405+00 
-.10745+00 
-.10915+00 
-.10915+00 
-.10755+00 
-.10455+00 
- .10005+00 
-.53975-01 
-.859 05-01 
-.75635-01 
-.62505-01 
-.45645-01 
-.26075-01 
-.43096-02 

,14325-01 


.10415-05 
-.93035-02 
15637-0 1 
-.20915-01 
-. 26075-C1 
-.31565-01 
-,3767»--01 
-.44555-01 
-.52245-01 
-.60655-01 
-.69555-01 
-. 786 O-T-0 1 
-.8735t>01 
-.95445-01 
-.10245+00 

-*ioeit+oo 

-.11245+00 
-.11515+00 
-.11655+00 
-.11655+00 
-.11555+00 
-.11335+00 
-.1097 r +Q0 
-.10425+00 
-.96255-01 
- 184885-0 1 
- • 694 4 1 " -0 1 
-.5-0365-01 
-.27535-01 
-.87935-02 


.58765-07 
-.91895-02 
-.15496-01 
*.20846-01 
-.26096-01 
-.31675-01 
-.37895-01 
-.44376-01 
-.52665-01 
-.61175-01 
-.7014 6-01 
-.79305-01 
-.88175-01 
-.96395-01 
-.10366+00 
-.10975+00 
-.11445+00 
-.11785+00 
-.12015+00 
-.12135+00 
-.12175+00 
-.12145+00 
-. 1200E+00 
-.11716+00 
-.11196+00 
-.10305+00 
-.89675-01 
-.72456-01 
-.50655-01 
-.32905-01 


-.2339E-0R 
-.91045-02 
-.15415-01 
-.20825-01 
-.26145-01 
-.31805-01 
-.38066-01 
-.45095-01 
-.52935-01 
-.61475-01 
-.70485-01 
-.79665-01 
-.88576-01 
-.9688E-01 
- • 1043E + 00 
-.11055+00 
-.11555+00 
-.11945+00 
-.12225+00 
-.12425+00 
-.12566+00 
-.12666+00 
-.12695+00 
-. 12606+00 
-.12335+00 
-.11726+00 
-.10615+00 
-.91515-01 
-.72335-01 
-.56035-01 


-.61975-09 
-.90646-02 
-.15365-01 
-.20775-01 
-.26085-01 
-.31735-01 
-.37995-01 
-.45016-01 
-.52845-01 
-.61385-01 
-.70406-01 
-.79605-01 
-.88566-01 
-.56965-01 
-.10456+00 
-.11095+00 
-.11625+00 
-.12045+00 
-.12365+00 
-.12616+00 
-.12835+00 
-.13026+00 
-. 1 317E + 0 0 
-.13245+00 
-. 1 31 5E + 0 0 
-.12805+00 
-.11.926 + 00 
-.10735+00 
-.91506-01 
-• 7785E-0 1 


- • 2739E-1 0 
-.9055E-02 
-.15335-01 
- .2069E-0 1 
-.25925-01 
-.3147E-01 
-.3763E-01 
-.4455E-01 
- . 523 2E-0 1 
-.6081E-01 
-.6982E-01 

- .7903E-01 
-.eeiiE-oi 
-.9661E-0 1 
-.10435+00 
-.11105+00 
-.11655+00 
-•1210E+00 
-.12465+00 
-.1276E+00 
-.13035+00 

'-.1329E+00 

- . 1354E + 0 0 
-.1373E+00 
-.1383E+00 
-.1362E+00 
-.1292 E+00 
-.1197E+0G 
-.1071 E+00 

- .9652E-0 1 



I z= - 

7 = 

1 .10005+01 • 

2 .20005+01 • 

3 .30005+01 • 

• 4 .4000E+01 < 

5 .50005+01 • 

6 .60005+01 ■ 

V .70005+01 ■ 

8 .60005+01 - 

9 .50005+01 - 

10 .1000E+02 - 

11 .11005+02 - 

12 .1200E+02 - 

13 • 1300E+02 - 

14 .14005+02 - 

15 . 1500E+02 - 

16 .16005+02 - 

17 • 17005+02 - 

18 . . 1P0.OE + 02 - 

19 . 15005+02 - 

20 .20005+02 - 

21 .21005+02 - 

22 .22005+02 - 

23 .23005+02 - 

24 .2400E+02 - 

25 .25005+02 - 

26 .2600E+02 - 

27 • 2700E*02 - 

28 .28005+02 - 

29 .25005+02 - 

30 >30005+02 - 


31 

•3100 £+02 
-.10235-10 
-.68025-02 
-.14515-01 
-.20185-01 
-.25405-01 
-.3096E-0 1 

• .37155-0 1 
•.44135-01 
-.51 93E-D 1 
-.60465-01 
-.69555-01 

78835-01 

-.88025-01 

• .96665 -0 1 
•.10455+00 
-.1113E+00 
-.11705 + 00 
-.12135+00 
.12565+00 
. 1289E + 0 0 
•13205+0 0 

• 13S1E+0 0 

• 1381E+0 0 

• 1 4D9E + 0 Q 

• 1430E + 0 0 
. 1421E+D 0 
. 13635+0 0 
. 1 2675+0 0 

• 1 187E+0 0 
.1112E+C 0 


.32 

. 3200E+02 
.25345-11 
-.57865-02 
-.14865+01 
-.20155-01 
-.25355-01 
-.30895-01 
-.37085-01 
-.4405E-01 
-.51845-01 
-.60375-01 
-.65465-01 
-.78775-01 
-.87995-01 
-.56665-01 
-.10465+00 
- . 1 1 15E + 0 0 
-.11735+00 
-.12225+00 
-.12625+00 
-.12575+00 
-.13315+00 
-.13655+00 
-.14005+00 
-.14345+00 
-.14635+00 
-.14605+00 
• . 1 4065 + 00 
•.134,15 + 00 
•*12605+00 
■• 121 “ 2£+00 


33 

.33005+02 
.10925-03 
-.86755-02 
-.1477E-01- 
-.20035-01 
- « 2522E + 0 1 ■ 
-.3076E-01 
-.36945-01 
-.43915-01 
-.51695-01 
-.60215-01 
- . 6 5305—0 1 

7065E-01 
-.87865-01 
-.96595-01 
-» 1 0455+00 
-.11155+00 
-.11745+00 
-.122 3J5 + 00 

- .12645 + 00 
-.13005+00 

- • 1 3355 + 0 0 
-• 1371E+00 
-• 1408E+00 
*• 1446E+00 
-.14785+00 
-. 1480E + 00 
-.14295+00 
■» 1371E+00 
• • 1 300E + 0 0 
••12785+00 


ST 57I(-li 


W-VEL 


7 = 

IY Y 

1 .1000G+01 

2 .20005+01 

3 .30005.01 

4 .46005*01 

5 .50005+01 

6 .60005+01 

7 .7000E+01 

8 • 8000E+01 

9 .50005*01 

10 • 1 C006 +02 

11 .11005+02 

12 . 1 2005+02 

13 .131)05*0? 

14 .14005-02 

15 .15005+02 

16 .160CL+C2 

17 .17005+02 

18 .16005+02 

19 .150GG+C2 

20 .20005+02 

21 .21005+02 

22 .22005+02 

23 .23005+02 

24 .2*005+02 

25 .25005+02 

26 .260 05+02 

27 .276)05+02 

26 .21005+02 

29 .29005+02 

30 .30005-02 


.10005+01 .20005+01 .20005+01 .40003+01 .50005+01 .60005.01 .70005+01 .80005+01 


.50005+01 


10 

.10005+02 


-.15625-13 
.23785-11 
-.48003-06 
-«6452t -08 
”•78261"- Is 
“•H«.i67£“Db 
-.9667C-Cb 
-•1 0261-07 
•1067E-07 
.10921 -07 
.11031-07 
•lOVvt -07 
.1G821-C7 

• lC52f>07 
•10101-07 
•99621-08 
•89361-Gb 
.6?<*6C-Di> 
•75*21 -Ob 
•68752-08 

• 6 3 0 9 - 0 b 
• 591 11-OH 
.57301-08 
.57811-Gb 
*5^94 L-CH 
.61561-08 

• 5 7 << r I ! - £ •. « 
.4C*cl-0 

• 4* 0G-..-1 1 
.54691-1 7 


.20481-03 

•64431-02 

• c 30 01-03 
.1 0351-02 
.96861-03 

• F 6771-03 
.14461-03 

-.601 11-02 
-.15001-02 
-.26981-02 
3V641-02 
-.13541-02 
- • 6 tt 0c-02 
-.81121-02 
-.9 359.-.-Q2 
104 4 ,->01 

-. 11321 -ai 

-.119F.1-0] 

-.i:7h£.-01 
-•12961-01 
-.13121-01 
-.13251-01 
-♦153&L-01 
1 35 51-01 
i 37 br.- oj 
-*1444>-01 
-.15761-01 
1 7 351-01 
i c .-«6t>ei 


• 4 1321-03 

• 1 2401-.02 
.17971-02 
.20911-02 
•20021-02 
.15271-02 

♦ 64436.-03 
•63841-03 
.23131-02 
. 4 3£S.->02 
.66151-02 

* 9 07 OF. - 02 
.1 llfil-Cl 

♦ 1 4061-01 
•lfc341-01 
.16331-01 

• 19971-01 

♦ 2124t>01 
.22171-01 

♦ 22b 4 - 0 1 
.23311-01 
.23681-01 
.24041-01 
.24431-01 
.2469F -Cl 

• 25601-01 
♦2716L-G1 

• 3009001 
.536*5 £-01 


.6511003 

♦18732-02 

•27081-02 

.313900? 

♦30432-02 

♦ 23971-02 

• 11621-02 
-• 39 4 .30 03 
-.2915002 

• 5711.002 

• c- 60200? 
.1234001 
. I5»7i -01 
•1931 r -0 1 
.22511-01 
♦ 2531001 

• 27 630G1 
•23441-01 
. 3079 5* - 01 
♦3177£-01 
.3249001 
.33061-01 
•3366c -01 
.34311-01 
•3511 O 03 
. 363 2;..- 0 1 

• 38 70.-C+1 

• 430SOC1 
*4846 ‘OOl 
.57271-01 


.8504003 

• 2400002 

• 3476E>02 ' 
.40701-02 

• 401 O c ‘- 02 

• 32 4 0 E- 02 
.17741-02 

-.40801-03 

• 32 / 3 h — 02 

• 6 7 4 9 F - 0 2 
•1071 l -Cl 
. i^ccr-oi 

. 19411-01 

.23731-01 

•2776E-01 

• 312 3001 
.3*23001 

• 365 3E- 01 
*3827001 
.3952001 
*4047001 
.412*001 
•4208E-01 

• 4295001 

• 4 4 0 4 F - 0 1 
. 4 5 6 r - C 1 
.4876001 
♦54201-01 
*60931-01 

• 7 2 2 6 F - 0 1 


.10601-02 

.29091-02 

.42361-02 

.50341-02 

*50541-02 

.42691-02 

• 2656E-C 2 
•20751-03 
.30701-02 
♦706 c ”-02 

• 1 1 6 6 0*0 1 
.1668001 
.21861-01 
.26961-01 
.31741-01 
.35971-01 

• 39511-0 1 
.42321-01 
.44451-01 
.46041-01 
.47261-01 
.48281-01 
.49331-01 
.50501-01 
•51821-01 
.5357- -01 
.56911-01 
.6311001 
.70231-01 

• 8.24 81-G1 


•12861-02 
.34011-02 
• .49811-02 
.60301-02 
♦62021-02* 
.54731-02 
•38091-02 
•11871-02 
-.23661-02 
-.67621-02 
-.11681-01 
-.17511-01 
-.23371-01 
-.29171-01 
-.34681-01 
3959 1- 0 t 
-•43761-01 
-.4711 1-01 
- • 4^631-01 
-.51661-01 
-.53151-01 
-.5444 1-01 
-.55651-01 
-.57021-01 
-.58441-01 
-.60041-01 
-.63411-01 
-.69691-01 
-.7674 £-01 
-.87561-01 


.15141-02 

.38991-02 

.57591-0? 

.70831-02 

.75011-02 

.69401-02 

..53451-02 

.26841-0? 

-.10311-0? 

-.57231-02 
.1126T-01 
.17441-01 
.23951-01 
.30531-01 
.36801-01 
.42491-01 
.47381-01 
.51581-01 
.54431-01 
.56801-01 
.'58551-01 
. 59851 -Cl 
.61001-01 
.62011-01 
• 6.3151-01 
.64211-01 
.67161-01 
•732^1-0! 
. 82 16 E - Cl 
.93911-01 


.17771-02 

• 45.921-02 
.69161-02 
.87041-02 
.94561-02 
• 9127E-G2 
•76571-02 
.49801-02 

• 1094f>02 
-♦ 3 C 50E-0 2 
-.10021-01 
-• 1 688Er 0 1 
-.24271-01 
-.31751-01 
-.39111-01 
-.45321-01 
--.51651-01 
-.56421-01 
-.60051-01 
-.62551-01 
-.64131-01 
- . 6485E-0 1 
-* 651 2E-0 1 
-.65011-01 
-.65151-01 
-.65671-01 
-.67931-01 
-.72551-01 
-.78781-01 
-.88701-01 


.20471-02 
• 6099E-02 
•9374E-02 
.11911-01 
.13291-01 
.13401-01 
.12151-01 
• 9 46 0 E-0 2 
.53121-02 
-.25101-03 
-.71191-02 
-.15041-01 
-.23701-01 
-.32671-01 
-.41511-01 
-.49561-01 
-.56541-01 
-.62031-01 
-.65781-01 
-.67861-01 
- .683 1 E-0 1 
-.67481-01 
6573E-01 
-• 6389 E-0 1 
- • 6 276E-0 1 
-.63431-01 
-.65301-01 
-.67991-01 
-.71971-01 
- • 7925E-0 1 




12 = 

11 

12 

13 

14 


Z = 

• 1 1 0 0E + 02 

* 1200E+02 

•1300E+02 

.1400E+02 

IY 

1 

Y 

. • lDO’OE + 0 1 

•2027E-D2 

• 2 225E-02 

\ 

. *2225 K -02 

.16336-02 

2 

•2000E+01 

.7390E-02 

• 1 078 E-0 1 

♦ 1 33 4 E-C 1 

.13526-01 

3 

.3000E>01 

• 1144E-ril 

•1649E-01 

•2045E-01 

.21236-01 

4 

•40006+01 

• 1467 1* - 01 

•2027Z-01 

* 2 472C-0 1 

.25896-01 

5 

•50006+01 

• 1 639E-0 1 

•21V9E-01 

* 2 64 4 C "0 1 

.27796-01 

6 

• 6 OOGE + 0 1 

• 1665E-D 1 

.2164c-01 

*2f c OE-01 

♦27536-01 

7 

• 7 000 6 + 01 

•1520E-01 

• 151liE“Ul 

.2316E-D1 

.25406-01 

8 

•8000E+01 

•12D2E-01 

•1461E-01 

• 184UE-01 

.21646-01 

9 

•90006+01 

•7020E-02 

•P037E-02 

« 1 17 5E-0 1 

•1648E-01 

• 10 

•.1 000 E + 02 

•2491E-03 

-.3741E-03 

•3535E-02 

•1024E-01 

11 

•1100E+02 

8139E-0 2 

-♦ 1 031E-0 1 

-.5787E-02 

.32986-02 

12 

• 1200 E+0 2 

-• 1784E-D 1 

-.2121E-01 

-•1E77E-01 

-.37776-02 

13 

• 1 300 6 + 02 

-.2832E-01 

-.3255E-01 

-.256GE-01 

-.10596-01 

14 

•14006+02 

-* • 39 1 1 l-0 1 

-•4340E-01 

-• 3462 E -01 

-.16306-01 

15 

• 1 500 E + 02 

-•4933E-01 

-•5305E-01 

-■4201E-01 

-•20586-01 

16 

« 1 6006+02 

- •5849L-IJ I 

-.6065E-C1- 

- *4710 c -01 

- • 2251 E-G 1 

17 

• 1 700 £ + 02 

— •6566E-13 1 

- • 6566E-01 

-.4^55 0L-01 

-.23096-01 

18 

• 1 800E+02 

-*70 48E-C 1 

-.f 775E-01 

-•4920 E -01 

-.21196-01 

19 

. 1900E+02 

-••7271E-C 1 

-•6713L-01 

-•4647E-01 

-.17676-01 

20 

• 2 000E + 02 

-.7254E-01 

-.f 413E-01 

-•41 9 IE -01 

-.13026-01 

21 

• 2 100E + 02 

7025E-G 1 

— ♦ E > 3 2 c. — 0 1 

- « 3620E-01 

-.78176-02 

22 

• 2 20 0 E+ 0 2 

- • 6648E-JD 1 

-.f 371E-01 

- • 3 U 2 4 fc. *- 0 1 

-.28516-02 

23 

* 2300E + 02 

-•6219E-01 

- • 4 8 2 3 E - 0 1 

- • 2 525E-0 1 

.85366-03 

24 

• • 2 400E.+ 02 

-•5834E-01 

-♦4409E-01 

-•2222E-01 

.22086-02 

25 

• 2 500E + 02 

-•5642E-01 

-• A 28’6E-01 

- • 2 263E-0 1 

.22926-03 

26 

•2800E+02 

-.5770E-C1 

- • 4*529E-Q 1 

- • 2 623L-C 1 

-.42736-02 

27 

• 2700E+02 

59J.3E-0 1 

-•4693E-01 

-.2686E-01 

-.62536-02 

sil 

• 2600E+02 

-.6009E-J31 

- • 4 766E-0 1 

- • 3 103E-0 1 

-.12416-01 

•2900E+02 

- •6285E-0 1 

-•5115E-01 

-•3587E-01 

-.19266-01 

30 

•3000E+02 

-•6929E-01 

-.f 917E-01 

-*45616-01 

-.29176-01 


J2 = 

21 

22 

23 

24 


7 . = 

•2100E+02 

•2200E+02 

•2300E+02 

.24006+02 

IY 

1 

Y 

• 1 OOOE+Ol 

•5471E-04 

if 779E-03 

.81076-03 

.69126-03 

2 

*2C00E-»'01 

• 224 4 E -0 2 

♦ 1 028E-02 

- • 16956-04 

-.49446-03 

3 

•3000 E +01 

* 5344 E-0 2 

•3 180 E "02 

.11026-02 

.46096-03 

4 

•4 OOOE+Ol 

• 1 0 94 E-0 1 

• f 22 IE -0 2 

.49896-02 

.40176-02 

5 

•5 OOOE+Ol 

•1915E-01 

• 1 60 2 E-0 1 

; 1 1246-01 

.98296-02 

6 

•6000E+01 

032E-0 1 

• 2 663E-0 1 

.20206-01 

.17946-01 

7 

^•7 OOOE+Ol 

• 4486E-0 1 

• 4 C83E-0 1 

• 3 192E-0 1 

.28406-01 

8 

*8 OOOE+Ol 

»6230E-0 1 

•f 805E-01 

•4647E-01 

•41146-01 

9 

•9 OOOE+Ol 

• 8302E-D 1 

. 7829E-01 

• 6 367E- 0 1 

.55966-01 

10 

•1000E+02 

• 1 063 E +0 0 

.1 010E+00 

•8310E-01 

.72446-01 

11 

• 1 1 0 0 E +.02 

• 1313E + G 0 

•1253E+00 

.10416+00 

.89996-01 

12 

•1200E+02 

• 1567E + 0 0 

..1500L+00 

.12576+00 

.10786+00 

13 

» 1300E+D2 

•1809E+00 

.1 737E+00 

. 14686+00 

.12516+00 

14 

•1 400E+02 

•2024E+00 

•1950E+0U 

.16646+00 

.14106+00 
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